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PREFACE 


The  Rand  Corporation  recently  developed  new  analytic  methods  for 
studying  the  transient  behavior  of  component-repair/ inventory  systems 
under  time-dependent  operational  demands  and  logistics  decisions  like 
those  that  might  be  experienced  in  wartime  (Hillestad  and  Carrillo, 

1980).  These  methods  have  culminated  in  the  Dyna-METRIC  mathematical 
model  described  in  this  report. 

Dyna-METRIC  evolved  through  a  series  of  projects  that  addressed 
policies  for  improving  aircraft  readiness  and  supportability . 

Initially,  the  dynamic  queueing  equations  were  developed  and  applied  to 
the  problem  of  determining  appropriate  levels  of  spare  engines  for  the 
C-141  aircraft  to  support  planned  wartime  surges  in  flying  activity 
(Berman,  Lippiatt,  and  Sims,  1978).  Later,  indentured  component 
features  and  aircraft  availability  performance  measures  were  added  to 
alternative  logistics  policies  for  carrier-based  aircraft 
squadrons  in  the  Defense  Resources  Management  Study  (Rice,  1979). 

Further  features  were  added  to  enable  the  study  of  multiple  bases  of 
aircraft  and  their  dependence  on  transportation  of  spare  parts,  in 
support  of  the  Project  AIR  FORCE  project  "Responsive  Intra-Theater 
Transportation  System  for  Spare  Parts."  The  model  was  then  extended  to 
include  test  equipment,  limited  service  capacity,  and  test-equipment 
failure  probability  in  the  Project  AIR  FORCE  project  "Supporting  Modern 
Tactical  Avionics."  Centralized  repair  and  resupply  capability  features 
were  added  to  allow  the  Dyna-METRIC  model  to  be  embedded  in  the 
experimental  Combat  Support  Capability  Management  System  being  tested  in 


the  Pacific  Air  Forces  (PACAF)  under  the  Project  AIR  FORCE  project 
"Combat  Support  Capability  Management  System."  Most  important,  the 
models  will  be  used  in  the  Air  Force's  Worldwide  Combat  Supplies 
Management  System  (CSMS).[1] 

Other  initial  implementations  of  the  Dyna-METRIC  model  have  been 
used  by  the  Ogden  Air  Logistics  Center,  Air  Force  Logistics  Command 
(AFLC),  for  studying  USAF  F-4  and  F-16  aircraft  readiness  and 
supportability ;  by  Headquarters  AFLC  for  F100  engine  evaluations;  and  by 
the  Tactical  Air  Command  to  study  the  effect  of  several  repair  and 
supply  strategies  on  the  readiness  and  deployability  of  F-15  tactical 
squadrons . 

This  report  is  intended  for  users  of  the  Dyna-METRIC  computer  model 
and  for  others  who  wish  to  understand  the  model's  underlying 
mathematics  Complex  theory  has  deliberately  been  avoided  to  facilitate 
this  understanding.  Hillestad  and  Carrillo  (1980)  describe  some  of  the 
mathemat."s  with  greater  theoretical  rigor.  The  report  does,  however, 
assume  in  the  reader  a  basic  knowledge  of  calculus  and  probability 
theory. 

The  report  describes  the  mathematical  approaches  to  modeling  the 
effects  of  spare  parts  supply,  component  repai  and  related  processes 
on  combat  capability.  It  does  not  describe  the  imp 1 ement at ion  of  any 
specific  version  of  Dyna-METRIC.  The  actual  implementation  is  described 
in  a  forthcoming  program  description  document  (Pyles  et  al.,  1982). 

[1]  Actually  there  are  several  programmed  versions  of  Dyna-METRIC. 
The  version  being  embedded  in  CSMS  is  referred  to  as  Dyna-METRIC,  while 
all  other  versions  are  called  RAMS  (Rand  Analytic  Models  of  Support). 

In  this  report,  the  author  refers  t*-*  the  entire  generic  class  of  models 
as  Dyna-METRIC. 
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The  publication  of  the  report  has  been  supported  by  the  "Combat 
Support  Capability  Management  System"  study  conducted  within  the  Project 
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SUMMARY 


This  report  describes  the  Dyna-METRIC  mathematical  model  for 
relating  aircraft  spare-parts  supply  levels  and  maintenance  capability 
to  material  readiness  of  aircraft.  A  key  characteristic  of  the  model  is 
its  ability  to  deal  with  the  dynamic  or  transient  demands  placed  on 
component  repair  and  inventory  support  caused  by  time  variables  in  a 
scenario  that  includes  sortie  rates,  mission  changes,  phased  arrival  of 
component  repair  resources,  interruptions  of  transportation,  and  the 
like. 

This  has  been  accomplished  by  implementing  a  set  of  analytic 
mathematical  equations  describing  the  dynamic  behavior  of  the  component 
repair  queueing  system.  The  set  also  includes  mathematical  models  of 
components  and  subcomponents  (indentures)  and  multiple  echelons  of 
repair  capability.  From  these  equations  the  model  can  compute  time- 
dependent  inventory  levels  (spare  parts  requirements)  and  assess  time- 
dependent  mission  readiness  of  the  aircraft  supported. 

The  text  develops  and  describes  the  mathematical  features  of  the 
model,  providing  insight  into  both  its  capabilities  and  the  inherent 


assumptions . 
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I.  INTRODUCTION 


Congress,  the  Department  of  Defense,  and  the  services  have  each 
recognized  the  importance  of  being  able  to  relate  logistics  resources 
and  policies  to  readiness .[ 1]  They  have  also  recognized  the  difficulty 
of  establishing  that  relationship.  For  example,  the  House  Armed 
Services  Committee  Readiness  Panel  reports  (1980)  that  "Current  DOD 
readiness  data  and  reports  are  not  particularly  helpful'  to  relating 
resources  to  future  readiness  ..."  The  DOD  Material  Readiness  Report 
(GAO,  1980)  states, 

The  Department  of  Defense  (DOD)  spends  billions  each  year  to 
maintain  the  readiness  of  its  weapon  systems  but  cannot 
accurately  project  how  much  readiness  a  dollar  will  buy  or 
determine  how  much  readiness  is  needed.  ...  To  date,  DOD 
has  made  little  progress  in  linking  funding  and  material 
readiness  and  has  not  achieved  an  adequate  material  readiness 
report  for  the  Congress.  Its  officials  told  GAO  that  the  size 
and  complexity  of  the  problem  has  seriously  hampered  their 
attempts . 

In  describing  the  difficulty  of  measuring  readiness,  Charles  W. 
Groover,  Deputy  Assistant  Secretary  of  Defense  (Program  Integration) 
MRA&L,  states  that  a  major  difficulty  in  managing  readiness  is,  "the 
absence  of  analytical  links  between  resources  and  readiness.  We  believe 
we  have  a  pretty  fair  understanding  of  how  the  logistics  system  operates 
to  support  our  combat  weapon  systems  and  equipment.  .  .  .  However,  the 
specific  functional  relationships  between  resources  applied  and  material 


[1]  Readiness  is  defined  as  the  state  of  availability  and 
preparedness  of  resources  with  respect  to  their  planned  wartime  mission. 
It  differs  from  capability,  a  measure  of  the  comparative  effectiveness 
of  different  resources  for  accomplishing  the  same  mission,  and  from 
effectiveness ,  which  relates  capability  to  its  application  against  enemy 
forces . 
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readiness  resulting  is  incredibly  complicated."  (Groover,  1977)  The 
Air  Force,  in  its  Logistics  Long  Range  Planning  Guide  (USAF,  1981), 
states  that  a  major  objective  is  to  "develop  a  means  to  better  identify 
and  assess  logistics  requirements  and  capability,  especially  as  these 
relate  to  execution  of  U.S.  contingency  plans." 

Currently,  readiness  reporting  by  the  services  consists  of  rating 
the  readiness  of  the  separate  groups  of  resources  necessary  to  support  a 
combat  unit:  supply,  fuel,  munitions,  aircraft,  operational  personnel, 
maintenance  personnel,  etc.  The  group  with  the  lowest  rating  defines 
the  overall  unit  readiness.  Unfortunately,  the  measures  used  for  the 
different  resources  are  not  comparable  and  do  not  equally  reflect  the 
ability  of  the  organization  to  perform  its  combat  mission.  For  example, 
the  supply  system  is  measured  by  its  ability  to  satisfy  demands  for 
spare  parts,  while  pilots  are  measured  by  the  number  who  are  proficient 
and  mission-ready.  To  achieve  a  balanced  state  of  readiness,  the 
services  need  techniques  for  describing  readiness  that  use  consistent 
measures  and  consider  the  interactions  among  resources. 

The  Dyna-METRIC  model  was  developed  to  study  and  predict  the 
readiness  of  groups  of  aircraft  squadrons  as  determined  by  a  major 
subset  of  logistics  resources,  namely,  those  associated  with  component 
repair  and  resupply.  Thus,  it  attempts  to  solve  two  of  the  problems 
mentioned  above:  combining  the  influence  of  several  types  of  support 
resources,  and  measuring  their  direct  effect  on  mission  readiness. 

Figure  1  illustrates  the  flow  of  aircraft  components  from  the  aircraft 
to  various  repair  facilities  to  on-hand  inventory,  and  from  this 
inventory  to  the  aircraft.  The  net  objective  of  this  system  is  to  avoid 
the  loss  of  aircraft  mission  capability  due  to  shortages  of  correctly 
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functioning  components  on  the  aircraft.  Clearly,  this  objective  can  be 
met  only  if  the  local  supply  of  those  components  exceeds  the  number  of 
components  tied  up  in  various  phases  of  repair  or  shipment.  A  snapshot 
in  time  would  show  some  components  awaiting  repair,  some  being  repaired, 
some  in  operable  condition  in  inventory,  some  on  the  way  to  and  from 
another  echelon  of  repair,  and  some  partially  repaired  but  awaiting 
spare  parts.  Each  of  these  states  is  a  "pipeline"  that  contains  some  of 
the  total  inventory  of  components .  When  there  are  not  enough  spare 
components  to  cover  each  of  these  pipelines,  "holes"  will  appear  in 
aircraft;  these  "holes"  may  or  may  not  affect  the  ability  of  the 
aircraft  to  perform  a  mission,  depending  on  the  mission  essentiality  of 
the  missing  component.  The  set  of  resources  considered  in  Dyna-METRIC 
includes  the  supply  of  spare  aircraft  components  at  bases  and  higher 
echelons,  transportation  resources,  and  the  personnel  and  equipment  used 
to  perform  repairs.  The  model  links  these  resources  to  the  unit's 
sortie  capability.  (Actually,  transportation  is  represented  by  delays 
and  interruptions,  and  personnel  and  equipment  are  represented  by  repair 
times  and  capacities  of  the  repair  process.) 

The  readiness  of  a  unit  to  perform  its  mission  depends  on  the 
availability  of  resources  to  support  a  highly  dynamic  flying  program  in 
the  face  of  delays  and  interruptions  of  logistics  support.  Current 
methods  for  calculating  required  amounts  of  resources  (such  as  spare 
parts  and  munitions)  use  "steady  state"  or  time-averaging  techniques 
that  do  not  account  explicitly  for  surges  in  sortie  demands  or 
variations  in  logistics  support.  Muckstadt  (1980)  has  shown  that  spare- 
parts  levels  determined  with  such  techniques  can  seriously  understate 
the  requirements  during  peak  periods  of  activity.  Berman,  Lippiatt,  and 


Fig.  1  --  Component  repair  and  the  readiness  assessment  problem 


Sims  (1978)  have  shown  that  logistics  strategies  (prepositioning  spare 
engines,  etc.)  that  are  not  apparent  in  the  steady-state  methods  become 
more  obvious  when  time  dependence  is  considered  explicitly.  A  key 
characteristic  of  the  Dyna-METRIC  model  is  its  ability  to  deal  directly 
with  the  transient  demands  placed  on  component  repair  and  inventory 
support  caused  by  dynamic  parameters  in  a  scenario  (sortie  rates, 
mission  changes,  phased  arrival  of  component  repair  resources, 
interruptions  of  transportation,  etc.).  It  does  this  by  implementing  a 
set  of  analytical  mathematical  equations  describing  the  dynamic  behavior 
of  the  component  repair  queueing  systems,  hence  the  term  "Dyna"  in  the 


title  of  the  model.  It  also  includes  mathematical  models  of  components 
and  subcomponents  (indentures)  and  multiple  echelons  of  repair 
capability.  From  these  equations,  it  can  compute  time-dependent 
inventory  levels  (spare  parts  requirements)  and  assess  the  mission 
readiness  of  the  aircraft  supported.  For  these  reasons,  we  borrowed  the 
name  METRIC  from  Sherbrooke  (1968);  this  is  an  acronym  for  Multi-Echelon 
Technique  for  Recoverable  Item  Control,  a  mathematical  model  used  to 
compute  optimal  inventory  requirements  for  steady-state  activity  levels. 

Hillestad  and  Carrillo  (1980)  give  the  theoretical  development  of 
the  dynamic  queueing  equations  that  form  the  heart  of  Dyna-METRIC,  as 
well  as  certain  performance  calculations  related  to  component  repair. 
This  report  reviews  those  results  and  describes  how  multiple  echelons, 
multiple  indentures,  and  requirements  optimization  can  be  incorporated 
in  a  dynamic  model  of  logistics  support.  In  describing  the  various 
aspects  of  the  model  below,  we  frequently  state  the  classical  steady- 
state  approach  prior  to  giving  the  time-dependent  result;  we  do  so  to 
illustrate  differences  in  the  approaches.  The  implementation  of 
Dyna-METRIC  is  described  by  Pyles  et  al.  (1982). 

Section  II  reviews  the  time-dependent  pipeline  equations.  Section 
III  describes  time -dependent  stockage  and  component-repair  measures  of 
performance.  Section  IV  combines  these  measures  to  give  aircraft 
capability  measures.  Section  V  introduces  the  pipeline  model  for 
indentured  components,  while  Section  VI  describes  the  pipeline  equations 
for  the  time-dependent,  multiple-echelon  model.  Section  VII  describes 
the  optimization  techniques  for  supply  requirements,  and  Section  VIII 
describes  the  approach  for  limited  service  facilities. 
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The  reader  should  be  aware  that  Dyna-METRIC  embodies  many  inherent 
assumptions  regarding  probability  distributions,  independence  of 
variables,  failure  and  repair  of  components,  and  allocation  of 
shortages.  Rather  than  list  them  out  of  context  and  without  evaluative 
statements,  we  discuss  the  assumptions  in  the  development  of  the  model 
in  the  following  sections. 
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II.  TIME -DEPENDENT  PIPELINES  AND  PROBABILITY  DISTRIBUTIONS 


Classical  steady-state  inventory  theory  provides  a  model  that 
describes  how  many  components  will  be  in  the  various  pipelines  of  a 
component -repair/ inventory  process  when  the  component  failure  rates  are 
driven  by  a  probability  distribution  that  is  not  time -dependent  (that 
is,  one  that  is  stationary).  Consider  the  case  of  a  single  type  of  com 
ponent  with  an  average  daily  failure  rate,  3,  and  an  average  repair 
time,  T.  Steady-state  theory  (Palm,  1938),  shows  that,  under  the 
assumptions  that  the  probability  distribution  of  repair  time  is  indepen¬ 
dent  of  the  failure  process,  and  that  ample  service  capacity  exists, [1] 
the  average  number  of  components  in  the  repair  pipeline  will  be[2] 


T 


Generally,  aircraft  components  are  assumed  to  fail 
the  number  of  flying  hours  so  that  we  can  break  d, 
demand  rate,  into  its  constituents: 


at  a  rate  based  on 
the  average  daily 


d  =  (average  failures  per  flying  hour)  x  (flying  hours/ 
V  sortie)  x  (average  number  of  sorties  per  day  per 

aircraft)  x  (number  of  aircraft)  x  (quantity  of  the 
components  on  the  aircraft)  x  (percentage  of  air¬ 
craft  with  the  component)  . 


[1]  The  assumption  of  ample  service  capacity  implies  that  no  com¬ 
ponent  awaits  service  and  prevents  the  steady-state  queue  from  becoming 
unbounded.  The  effect  of  limited  service  capacity  is  discussed  further 
in  Sec.  VIII. 

aware  that  we  have  changed  the  usual  mean 
which  states  the  same  result  as  L  =  XW. 


[2] 

ing  of  X 
Thus  L  = 


The  reader  should  be 
in  Little's  Theorem, 
*ss.  »  -  3.  «  =  T. 


i 


I 
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With  the  further  assumption  that  d  has  a  Poisson  probability  distribu¬ 
tion  the  average  number  of  components  in  the  steady-state  pipeline,  Xgs, 
actually  represents  the  mean  value  of  a  Poisson  distribution  as  well. 

By  Palm’s  Theorem,  the  probability  that  there  are  k  components  in  the 
pipeline  at  any  point  in  time  is  given  by 

.  k  -X 
X  ess 

P(k  in  pipeline)  =  -  ss-  ^ - 

The  Poisson  distribution  arises  when  the  mean  time  (in  flying  hours) 
between  failure  of  components  follows  an  exponential  probability  distri¬ 
bution,  and  is  frequently  borne  out  in  data  analysis  of  failure  rates. 
There  are  also  strong  theoretical  reasons  to  expect  Poisson  failures. [3] 
Time -dependent  flying  scenarios  cannot  be  considered  correctly  with 
this  approach.  It  is  necessary  to  use  the  average  number  of  aircraft 
and  average  number  of  sorties  per  day  per  aircraft  rather  than  the 
time-dependent  values.  In  addition,  time -dependent  repair  capability 
cannot  be  considered  because  the  average  repair  time,  T,  must  be  used  in 
the  steady-state  model. 

The  dynamic  model  used  in  Dyna-METRIC  relieves  these  assumptions. 

In  this  case  we  let  the  daily  demand  rate,  d(t),  be  a  function  of  time 
so  that 

d(t)  =  (failures  per  flying  hour)  *  (flying  hours/sortie  at 
time  t)  x  (number  of  sorties  per  day  per  aircraft  at 
time  t)  x  (number  of  aircraft  at  time  t)  x  (quantity 
of  the  component  on  the  aircraft;  x  (percentage  of 
aircraft  with  the  component)  . 

[3]  See,  e.g.,  Samuel  Karlin  and  H.  M.  Taylor,  A  First  Course  in 
Stochastic  Processes ,  2d  ed. ,  Academic  Press,  New  York,  1975,  pp.  221- 
228.  ~ . 
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This  expression  for  d(t)  includes  severe  1  variables  that  might  be 
expected  to  change.  For  example,  the  number  of  aircraft  can  change  in 
time  according  to  the  time  sequence  of  deployment  or  because  of  aircraft 
attrition.  The  number  of  sorties  per  day  per  aircraft  changes  as  a 
result  of  programmed  changes  in  flying  rates,  and  the  flying  hours  per 
sortie  change  as  missions  change. 

In  place  of  a  constant  average  repair  time,  T,  the  dynamic  model 
uses  the  probability  that  a  component  entering  repair  at  time  s  is  still 
in  repair  at  time  t.  This  probability  function,  F(t,s),  is  called  the 
repair  function.  It  is  defined  by 

_  ( Component  entering  at  s  is 

F(t,s)  -  Prob  } still  in  repair  at  t 

Repair  time  >  t-s  when) 
started  at  s  f 


A  few  illustrations  will  show  that  this  function  is  relatively  sim¬ 
ple  to  obtain  for  most  uses,  but  also  has  considerable  power  to 
represent  certain  time  dependencies  in  component  repair.  (It  is  assumed 
to  be  independent  of  the  failure  process,  however.) 
i.  Constant  or  fixed  repair  time,  T. 

!1  if  t  -  8  <  T 
0  if  t  -  s  >  T 


(That  is,  if  the  time  between  component  arrival,  s,  for  repair 
and  the  current  time,  t,  is  less  than  the  constant  repair  time, 
T,  then  the  probability  that  it  will  stil]  be  in  repair  is  1; 
otherwise  it  is  0). 


(5) 
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ii.  Exponentially  distributed  repair  time  with  average  T. 

‘  t~8 

;  F(t,s)  -  e~  T  (6) 


iii.  No  repair  capability  until  time,  t,  with  exponentially  dis¬ 
tributed  repair  time  after  t. 


if  s  <  t  <  t 


if  T  <  S  <  t 


(7) 


iv.  Fixed  transportation  lag,  S,  with  exponentially  distributed 
repair  time  after  transportation  occurs. 


if  t  -  s  <  S 


F(t,s) 


t-(s+S) 

T 


if  t  -  8  >  S 


(8) 


v.  Exponentially  distributed  repair  time  which  changes  at 
time,  t,  (from  to  average  repair  time) 


F(t,s) 


t— 8 


t-T 


t-s 


T-S 


if  t  <  T 


if  8  <  T  <  t 


if  T  <  8  <  t 


(9) 
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Of  course,  many  other  repair  functions  could  be  modeled  provided  they 
are  independent  of  the  demand  function.  Section  VIII  shows  how  depen¬ 
dent  repair  functions  may  be  approximated.  We  will  now  describe  how 
Dyna-METRIC  combines  the  repair  and  demand  functions  to  determine  the 
average  number  of  parts  in  the  pipeline.  Consider  only  those  components 
that  arrived  in  an  interval  of  time,  As,  centered  at  time  s.  The  number 
expected  in  the  repair  pipeline  at  time  t  will  then  be  given  by 


AA(t,s)  *  d(s)  x  F(t,s)  x  As  , 


(10) 


where  A\(t,s) 

d(s) 

F(t,s) 

As 


expected  number  of  components  in  the  repair  pipeline 
at  time  t  that  arrived  during  the  interval  around  s; 

daily  failure  rate  at  time  s; 

probability  of  component  not  out  of  repair  by  time  t;  and 
interval  of  time  centered  at  s. 


If  we  assume  that  the  number  of  failures  arriving  in  the  interval  As  is 
independent  of  the  number  of  failures  arriving  in  similar  intervals  cen¬ 
tered  at  other  times  other  than  s[4]  and  that  the  repair  probability 
function  is  independent  of  the  probability  distribution  generating  the 
demand  rate,  we  can  sum  the  contributions  of  all  intervals  to  obtain 


[4]  This  is  called  the  independent  increment  assumption  and  is  true 
for  demands  governed  by  a  Poisson  process.  Actually,  there  are  reasons 
that  component  repair  demands  may  not  have  independent  increments.  If 
aircraft  are  decreasing  in  availability  because  of  component  failures 
and  there  is  a  resulting  decrease  in  sorties  flown,  then  the  later  sor¬ 
ties  (and  hence  component  failures)  are  a  function  of  previous  component 
removals.  This  will  cause  the  result  above  to  be  an  overstatement  of 
the  number  of  components  in  the  pipeline.  For  most  users  of  Dyna-METRIC 
this  will  not  seriously  affect  the  answer;  but  for  some  cases  with  seri¬ 
ous  shortages  of  supply  resources,  a  second  iteration  with  reduced  sor¬ 
ties  based  on  the  previous  iteration  will  give  a  more  accurate  answer. 
Generally,  the  overstatement  of  pipelines  that  may  occur  will  give  con¬ 
servative  answers  to  stockage  requirements  and  to  capability  with  a 
fixed  level  of  stock. 
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X(t)  -  Z  AA(t,8) 

8<t 

(ID 

=  d(s)  .  F(t,s)As 
8<t 

If  we  make  ds  arbitrarily  small,  we  arrive  at  the  integral  expression 
for  X(t). 


X(t)  =  J  d(s)F(t,s)ds  (12) 

0 

This  integral  represents  the  average  number  of  components  in  the  repair 
pipeline  at  time  t.  Hillestad  and  Carrillo  (1980)  show  that,  with  the 
additional  assumption  that  the  component  failure  probability  distribu¬ 
tion  is  Poisson,  X(t)  is  the  mean  of  a  nonhomogeneous  (time  varying) 
Poisson  process.  That  is,  the  probability  of  k  components  in  repair  at 
time  t  is 


P(k) 


X(t)k 


e"  *(t) 
k! 


(13) 


where 


X(t) 


d(s)F(s,t)ds  . 


(14) 


The  integration  from  0  to  t  assumes  no  flying  before  0  (or  at  least 
F(s,t)  =  0  for  s  <  0).  Superposition  allows  consideration  of  flying 
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during  this  period.  Peacetime  or  steady-state  pipeline  quantity  preced¬ 
ing  t  =  0  can  be  considered  by  computing  X (0)  •  F(t,0),  the  quantity 
of  the  pipeline  remaining  at  time  t.  Under  other  assumptions  it  is  pos¬ 
sible  for  X (t)  to  represent  the  mean  value  of  other  probability  distri¬ 
butions.  Compound  Poisson  distributions  also  satisfy  the  independent 
increment  assumption  and  lead  to  pipeline  distributions  such  as  the 
negative  binomial  distribution  given  by 


POO 


(r  +  k  -  1)!  p* 
(r  -  1)  ! k!  '  r+k 

q 


(15) 


where  q  =  variance  to  mean  ration  (q  >  1) 


r 


X(t) 
q  -  1 


(16) 


p  =  q  -  1 


(17) 


and  X(t)  is  given  by  the  equation  shown  above.  This  distribution  can 
arise  when  groups  of  failures  occur  at  instants  of  time  (e.g. ,  due  to 
the  arrival  of  groups  of  aircraft)  and  the  number  of  instants  of  failure 
(number  of  batches  of  failures)  in  a  given  time  interval  is  governed  by 
a  nonhomogeneous  Poisson  distribution.  When  the  distribution  of  the 
number  of  failures  in  any  group  (at  any  instant  of  time)  is  given  by  a 
logarithmic  distribution,  the  resulting  distribution  of  the  pipeline's 
quantity  is  negative  binomial.  (The  Poisson  distribution  allows  only 
one  failure  at  each  instant.)  The  additional  information  necessary  to 
use  this  distribution  is  q,  the  var iance-to-mean  ratio.  This  must  be 
determined  by  statistical  estimation  of  the  variance  of  historical 
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failure  data.  (The  Poisson  distribution  has  a  variance-to-mean  ratio  of 
1  and  does  not  need  this  variable.) 

An  important  aspect  of  estimating  pipeline  quantities  is  the  abil¬ 
ity  to  separate  pipelines  with  dissimilar  demand  rates  and  repair  func¬ 
tions  for  ease  of  calculation,  and  then  to  sum  the  result  to  obtain  the 
average  total  quantity  of  components  tied  up  in  all  of  the  pipelines. 
This  is  easiest  when  the  components  flowing  into  each  pipeline  arise 
from  a  common  source  distribution.  For  example,  suppose  that  demands 
for  component  repairs  arise  at  an  aircraft  squadron  based  on  the  flying 
activity  of  the  squadron,  and  that  there  is  a  probability  p  that  a 
failed  component  can  be  fixed  locally  and  1  -  p  that  it  must  be  shipped 
elsewhere  to  be  repaired.  Further,  assume  that  the  local  repair  has 
probability  function  F^(s,t)  and  that  remote  repair  has  probability 

function  F2(s,t)  (including  transportation  time).  Then  we  have  that 

x(t)  -  j"  d(s) [pF^(s,t)  +  (1  -  p)F2(s,t)]ds  (18) 

0 


or  that 


t 

d( s)pF^ (s , t)ds 


(19) 


*2<t) 


t 

d(s)(l  -  p)F2(s,t)ds 


(20) 


A(t)  -  X1(t)  +  A2(t)  • 


(21) 
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Thus,  it  is  possible  to  separate  and  recombine  the  average  values  for 
the  various  component  repair  pipelines  in  the  dynamic  model  just  as  in 
the  steady-state  model.  This  is  true  regardless  of  the  probability  dis 
tributions  of  different  parts  of  the  pipeline.  However,  when  the  dis¬ 
tributions  are  of  different  form  or  have  differing  variance-to-mean 
ratios,  the  resulting  total  distributions  can  be  complicated.  When  the 
distributions  in  each  pipeline  are  independent,  we  can  determine  the 
resulting  distribution  by  the  more  complex  process  of  convolution. 

This  section  will  conclude  witn  an  illustration  of  the  dynamic 
pipeline  calculation  and  time-dependent  probabilities.  Consider  an  air 
craft  squadron  at  a  single  location  with: 

N  (s)  =  aircraft  at  time  s, 

Si 

D(s)  =  sorties  per  aircraft  demanded  at  time  s,  and 
FH(s)  =  flying  hours  per  sortie  at  time  s. 

A  component  on  the  aircraft  has  the  following  characteristics: 

m  =  failures  rer  flying  hour; 
q  =  quantity  per  aircraft. 

The  demand  function  for  this  squadron  for  this  component  is  then: 

d(s)  =  N  (s)  «  D(s)  •  FH(s)  .  m  .  q 
a 

Each  failure  of  the  component  has  a  probability,  p,  of  being  repaired 
locally  and  1  -  p  of  being  repaired  remotely.  If  it  is  repaired 
locally,  then  it  can  be  repaired  in  an  exponentially  distributed  time 
with  mean  T^  except  for  the  first  days  during  which  no  repair  capa- 


(22) 
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billty  exists.  If  it  is  repaired  remotely,  an  order  is  placed  for  a  new 

component,  which  is  immediately  shipped  and  arrives  after  a  fixed  tran- 
o 

sportation  delay  T  .  The  number  of  aircraft  and  flying  hours  per  sortie 
will  be  assumed  constant  in  this  example  so  that  NA(s)  =  Na  and  FH(s)  = 
FH. 

The  demanded  flying  program  will  have  one  step  at  time  zero  to 
sorties  per  day,  and  a  second  step  to  a  lower  number  of  sorties  D2  per 

12  2 

day  at  (We  will  assume  that  ip  >  x^,  that  T  >  T  ,  and  that  T  > 

Xq.)  The  pipeline  integral  for  the  component  is  then 


A(t) 


d(s)[pF1(s,t)  +  (1  -  p)F2(s,t)]ds  , 


(23) 


where 


d(s) 


FH 


m 


FH 


m 


s  <  x. 


8  ^  td 


(24) 


where 


_1 

F  (s,t) 


t-s 


‘  <  tr 


8  <  tr  <  c 


(23) 


Tr  <  8  <  t 
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F2(s,t) 


if  t  -  s  <  T2 


if  t  -  s  >  T2 


(26) 


Let  XX(t)  be  the  average  local  repair  pipeline  quantity  at  time  t,  and 
2 

let  X  (t)  be  the  average  remote  repair  pipeline  at  time  t.  Then, 


xx(t) 


d(s)  pF^(s,t)ds 


(27) 


and 


x2(t) 


d(s)(l  -  p)F2(s,t)ds 


(28) 


Performing  the  integration  we  obtain  the  average  pipeline  quantities 
given  by 


Average  number  of  components 
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Figure  2  illustrates  these  results. 


Fig.  2— Illustration  of  average  pipeline  calculations 


The  addition  of  the  probability  distribution  allows  us  to  determine 
the  probability  of  various  pipeline  quantities.  For  example,  if  the 
demand  distribution  is  Poisson,  the  pipeline  distribution  will  be  Pois¬ 
son,  as  described  earlier.  We  can  determine,  for  this  example,  the  pro¬ 
bability  that  no  more  than  K  components  will  be  in  the  local  repair 
pipeline  using  the  following  calculations: 

p(no.  in  local  pipeline  <  K)  ■ 


K 

E 

k-0 


X1(t)ke"Xl(t) 

k! 


(33) 


Quantity 


Note  that  there  are  two  tcalei 


Fig.  3— The  probability  that  no  more  than  K  components  are  in  local  repair 


Another  important  use  of  the  probability  distribution  is  to  deter¬ 
mine  confidence  intervals.  Suppose,  in  the  example,  we  wish  to  deter¬ 
mine  the  number  of  components  in  local  repair  that  we  are  quite  sure 
will  not  be  exceeded.  We  might  quantify  this  by  looking  for  the  number 
that  will  not  be  exceeded  with  90  percent  probability.  That  is,  we 
desire  90  percent  confidence  (or  .9  probability)  that  the  number  will 
not  be  exceeded. 

Let  L(t)  be  the  number  of  components  that  satisfies  this  confi¬ 
dence.  The  value  of  L(t)  can  be  determined  by  summing  the  Poisson  terms 
until  they  equal  or  exceed  .9  at  each  time  instant  of  interest.  Thus, 


Number  of  components 
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we  find  L(t)  by  solving  the  following  equation  for  the  smallest  value  of 
L(t): 

Xl(t)ke"Xl(t) 

.9  =  22  k! 

k=0 

This  is  illustrated  in  Fig.  4. 


Fig.  4— Illustration  of  confidence  level 
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III.  TIME -DEPENDENT  COMPONENT  PERFORMANCE  MEASURES 


Section  II  described  the  calculation  of  time-dependent  average  com¬ 
ponent  pipeline  quantities  and  their  associated  probability  distribu¬ 
tions.  Combining  these  pipeline  quantities  with  the  supply  levels  at 
the  same  instant  of  time  allows  the  determination  of  various  measures 
describing  the  availability  or  shortages  of  individual  components  at  the 
aircraft.  These  are  the  supply  or  component -repair  measures  used  by  the 
services  to  describe  supply  system  performance,  except  that  in  this  case 
time  dependence  is  included.  They  are  also  similar  to  the  measures  used 
by  the  services  to  determine  the  required  spare-parts  levels.  These 
measures  do  not  completely  describe  aircraft  mission  readiness,  since 
they  do  not  describe  the  combined  loss  of  capability  that  shortages  of 
groups  of  components  inflict  on  the  aircraft  when  cannibalization  (con¬ 
solidation  of  shortages)  and  mission  essentiality  are  considered.  They 
do  form  a  basis  for  higher- level  measures,  and  they  are  also  useful  in 
identifying  problem  components  once  overall  system  performance  has  been 
determined  to  be  inadequate.  The  next  section  will  deal  with  the  joint 
effect  of  the  components. 

The  component  performance  measures  are  determined  for  a  given  sup¬ 
ply  level,  S(t).  This  level  is  provided  to  protect  the  aircraft  squad¬ 
ron  from  shortages  due  to  components  in  repair  or  on  order.  When  the 
number  of  components  in  repair  plus  on  order  at  time  t  exceeds  S(t), 
then  the  system  is  said  to  be  in  a  "backorder"  state.  When  components 
are  backordered,  the  aircraft  will  be  missing  components,  possibly 
degrading  the  mission  or  sortie  capability. 


The  component  measures  typically  computed  by  the  Dyna-METRIC  model 


are : 


R(t)  5  Ready  rate  at  time  t--the  probability  that  an  item 
observed  at  time  t  has  no  backorders. 

FR(t)  =  Fill  rate  at  time  t--the  probability  that  a  demand 
at  time  t  can  be  filled  immediately  from  stock  on 
hand. 


EB(t)  =  Expected  back  orders--the  average  number  of 
shortages  of  a  component  at  time  t. 

VBO(t)  =  Variance  of  the  backorders,  a  measure  of  the 
random  variation  of  backorders. 


DT(t)  =  Average  cumulative  demands  by  time  t. 


The  last  measure  is  derived  from  the  daily  demand  rate,  d(s),  as: 


DT(t) 


d( s)ds 


(1) 


The  remaining  measures  use  the  average  pipeline  quantity,  X(t),  the 
stock  level,  S(t),  and  the  probability  distribution  P(k/X(t))[l]  chosen 
from  the  Poisson  distribution  (var iance-to-mean  ratio  =  1.0)  or  the 
negative  binomial  (var iance-to-mean  ratio  >  1.0)  distribution. 

The  readj  rate  is  given  by 


S(t 

t)  =  ^  P(k/X(t) ) 


(2) 


[1]  The  notation  P(k/X(t))  is  introduced  here  to  indicate  the  time 
dependence  and  mean  of  the  probability  distribution.  It  is,  however, 
the  same  probability  as  described  by  P(k)  in  Sec.  II. 


-26- 


The  interpretation  of  this  expression  is  as  follows:  P(k/X(t))  is  the 
probability  of  there  being  exactly  k  components  in  the  pipeline  when  the 
average  number  is  X(t).  P(0/X(t))  is  the  probability  of  zero  in  the 

pipeline.  The  sum  of  P(0/X(t)),  P(l/\(t)),  up  to  P(S/X(t))  is  the 

probability  that  S  or  fewer  components  will  be  in  the  pipeline.  This  is 
the  probability  that  the  pipeline  will  contain  no  more  than  the  stock 
level  and  is  therefore  the  probability  that  there  will  be  no  backorders. 

The  fill  rate  is  given  by 


S(t)-1 

FR(t)  -  V  P(k/X(t))  .  (3) 

k-0 


The  interpretation  of  this  expression  is  as  follows:  The  fill  rate  is 
the  probability  that  a  component  will  be  available  when  a  demand  is 
placed.  It  is  therefore  the  probability  that  demands  have  left  at  least 
1  component  available  and  is  therefore  the  sum  of  the  probabilities  of 
demands  less  than  the  stock  level. 

Expected  backorders  are  given  by 


EB(t)  -  E  (k  -  S(t)  P(k/X(t)> 
k-S(t)+l 


SCt) 

-  A(t)  -  S(t)  +  £  (S(t)  -  k)  P(k/A(t)) 

k-0 


(4) 


The  interpretation  of  the  first  expression  is  as  follows:  There  are  no 
backorders  if  stock  equals  or  exceeds  the  demand.  Therefore,  for  k  less 
than  or  equal  to  S(t),  the  backorders  are  zero.  For  k  greater  than 
S(t),  the  backorders  are  merely  (k  -  S(t)).  The  probability  of  any 
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demand  level,  k.  is  P(k/X(t)),  and  the  expected  value  of  the  backorders 
is  merely  the  product  of  the  various  values  the  backorders  can  take  on 
times  the  probability  of  a  demand  at  that  given  value.  The  second 
expression  for  the  expected  value  is  more  useful  computationally;  it  is 
obtained  by  using  the  facts  that 


£  kP(k/A(t) )  -  X(t)  and  (5) 

k-0 


£  S(t)P(k/A(t))  -  S(t)  .  (6) 

k-0 


The  variance  in  backorders  is  given  by 


00 

VB( t)  £  [k  -  S(t)]2  P(k/A(t))  -  [EB( t) ] 2  .  (7) 

k-S(t)+l 


For  the  Poisson  Distribution  this  is  given  also  by  the  computationally 
tractable  expression, 

VB( t)  -  A(t)  +  [ A( t)  -  S(t) ] 2  -  [EB( t) ]2 


(8) 


[k  -  S(t) ] 2  P(k/A(t))  . 
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This  expression  was  derived  using  the  previous  results  and  the  fact 
that,  for  the  Poisson  distribution, 


X)(k  -  X(t))2  P(k/X(t))  =  X ( t)  .  (9) 

k=0 


These  measures  can  be  computed  for  each  component  for  each  location 
and  time  of  interest.  They  provide  a  detailed  picture  of  the 
component's  status.  Figure  5  illustrates  the  component  measures  for  a 
case  in  which  the  average  pipeline  quantity  undergoes  a  buildup  in  time. 


IV.  TIME -DEPENDENT  SYSTEM  PERFORMANCE  MEASURES 


Dyna-METRIC  also  forecasts  the  effect  of  component  shortages  (due 
to  shortage  of  spares,  inadequate  component-repair  capability,  etc.)  on 
the  number  of  mission-capable  aircraft  and  consequent  ability  to  gen¬ 
erate  missions.  These  measures  show  the  combined  effect  of  the  com¬ 
ponent  shortages  on  aircraft  availability,  both  with  and  without  "canni¬ 
balization"  to  minimize  the  number  of  aircraft  with  "holes."  Under  can¬ 
nibalization,  the  model  assumes  the  ability  to  instantly  consolidate 
shortages  onto  the  smallest  number  of  airframes.  In  actuality  this 
would  be  done  only  when  the  aircraft  are  needed.  Thus,  this  measure  may 
overstate  the  number  of  aircraft  that  might  be  available  in  unstressed 
conditions  (when  all  aircraft  are  not  needed).  On  the  other  hand,  the 
measures  for  conditions  of  no  cannibalization  may  understate  aircraft 
availability,  because  some  cannibalization  almost  always  occurs  in  prac¬ 
tice.  The  true  answer  probably  lies  somewhere  in  between  these  bounds, 
closer  to  the  full  cannibalization  values  in  stressful  conditions  and 
closer  to  the  no-cannibalization  numbers  under  less  stressful  condi¬ 
tions.  The  average  total  number  of  backorders  is  a  useful  system  meas¬ 
ure  that  describes  the  total  number  of  holes  in  aircraft. 

Mission  essentiality  of  components  can  be  accommodated  by  sets  of 
indicators  for  each  component  showing  whether  or  not  the  component  is 
required  for  a  given,  finite  set  of  missions.  Shortages  of  components 
that  are  not  considered  mission-essential  at  a  point  in  time  are  then 
not  included  in  determining  the  aircraft  performance  measures. 
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The  following  set  of  system  performance  measure  equations  is 
repeated  from  Hillestad  and  Carrillo  (1980): 

1.  Average  number  of  systems  NMC  (not  mission  capable)  without 
cannibalization,  EN(t). 

This  measure  gives  the  average  combined  effect  of  recoverable  item 
shortages  on  the  aircraft  that  those  items  support.  Given  that  NA(t) 
represents  the  number  of  aircraft  that  are  supported  at  time  t,  and  that 
there  are  N  types  of  recoverable  items  required  on  each  aircraft,  assume 
that  the  shortages  of  any  one  of  these  items  will  make  the  system  nono- 
perational.  Also,  assume' that  shortages  of  components  cannot  be  conso¬ 
lidated  among  the  systems. 

The  probability  that  an  arbitrary  aircraft  is  missing  component 
type  i  when  there  are  k  shortages  of  i  across  the  fleet  of  NA(t)  air¬ 
craft  is  f  1  ] 


k 

NA(  t) 


(1) 


and  therefore,  the  probability  that  an  arbitrary  aircraft  has  a  shortage 
of  item  i  is  given  by  summing  across  the  possible  values  of  k  =  £  -  s.(t) 
times  the  probability  they  occur: 

”  <*  -  S  (t))  EB  (t) 

St(t)+1  P(£Ai(t))  =  NA(t)  ’  (2) 

where  EBi(t)  is  the  expected  backorders  on  component  i.  Assuming  that 
failures  are  independent,  the  probability  that  an  aircraft  is  non- 
mission-capable  due  to  shortages  of  some  item  is 


[1]  This  assumes  that  there  is  at  most  1  of  each  component  of  type 
i  on  any  aircraft.  The  modification  for  multiple  occurrences  on  an  air¬ 
craft  will  be  discussed  shortly. 
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And,  finally,  the  expected  number  of  non-mission-capable  aircraft  at 
time  t  is 


EN(t)  =  NA( t)  1 


N  /  EB  (t)\ 

/A  v  '“w); 


This  derivation  required  the  assumption  that  only  one  component  of  a 
given  type  was  on  each  aircraft.  If  this  is  not  the  case,  let  be  the 
quantity  of  item  i  per  system.  If  no  consolidation  of  shortages  occurs, 


EN(t)  =  NA(t)  1  -  (I  Pr  (acft  have  no  shortages  of  i) 

i=l 


NA(t)  1  - 


N  NA(  t)  ♦( 

n  £ 

i=l  y=0 


Pr  (acft  have  no  shortages  of  i  when 
y  shortages  of  part  i  exist  in  the 
squadron) 

(5) 


=  NA(t)  1  - 


N  NA( t) «Q . 

II  £  1 

i=l  y=Q 


NT’i 

/Q  -NA(t)  \  PBi^y> 

(  ) 


where  PBi(y)  is  the  probability  that  item  i  has  y  shortages  at  time  t 
(the  probability  distribution  of  backorders  of  component  i). 


si(t) 

£  P(k/X  (t)) 
k-0 


y  *  0 


PB^y) 


P(k  +  S(t)/A1(t)) 


y  >  0 
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2.  Average  number  of  aircraft  NMC  (non -miss ion -capable)  with  an 

instantaneous  cannibalization  policy.  EN  (t). 

-  -  -  c 

In  this  measure  it  is  assumed  that  shortages  of  all  components  are 
consolidated  to  make  the  smallest  number  of  aircraft  non-mission- 
capable.  Let  P1(j)  be  the  probability  that  shortages  of  the  ith  item 
are  less  than  or  equal  to  j.  Then 


si(t)+j 

P^j)  =  £  P(k/X  (t))  .  (7) 

k=0 


Let  P(j)  be  the  probability  that  the  number  of  non-mission-capable  air¬ 
craft  is  less  than  or  equal  to  j.  Then 

N  i 

P(j)  -  II  P  (j)  •  (8) 

i=l 


If  there  is  more  than  one  item  of  a  type  on  each  aircraft,  we  again 
employ  as  the  number  of  item  i  per  aircraft  and  obtain 

N  i 

P(j)  =  II  P  (Q,  •  j)  •  (9) 

i=l 

This  measure  is  important  by  itself  as  well  as  being  a  step  toward 
obtaining  the  average  number  of  NMC  aircraft  unde/  full  cannibalization. 
It  is  used  as  a  performance  constraint  equation  in  the  optimization  pro¬ 
cedure  in  Sec.  VIII.  Ebeling  (1978)  refers  to  this  distribution  as  the 
"Operational  Ready  Rate,"  and  argues  that  it  is  the  most  operationally 


1 


t 
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oriented  performance  measure  among  those  commonly  considered  for  a 
spare-parts  objective  function. 

The  expected  number  of  NMC  aircraft  with  full  cannibalization  is 
then  (using  the  fact  that  the  expected  value  of  a  nonnegative  distribu¬ 
tion  is  equal  to  the  sum  of  the  complement  of  that  distribution): 

NAm-1 

EN  (t)  =  2l  t1  -  pU>]  •  (10) 

C  j=0 


« 


The  NMC  distribution  function  with  full  cannibalization  is 

PNjCt)  =  PC j)  "  P(j  -  1)  ,  (11) 

which  then  can  be  used  to  give  the  variance  in  the  number  of  NMC  air¬ 
craft  , 


VN(  t) 


NA(t)  , 

E  i  pN/t) 

J-l  J 


-  ENc(t) 


(12) 


3.  Total  expected  backorders . 

This  measure  is  the  sum  of  the  individual  component  backorders, 
BO^Ct).  Thus, 


N 

EB(t)  -  £  EB  (t)  . 
i-1 


(13) 
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i 


This  is  commonly  used  as  an  optimization  objective  or  constraint  in 
determining  supply  levels.  It  is  a  very  conservative  estimate  of  opera 
tional  performance,  as  can  be  seen  by  expanding  the  expected  NMC  expres 
sion  with  no  cannibalization  and  assuming  =  1  for  all  i. 


EN(t)  =  NA( t) 


1  - 


N  EB.(t) 

E  1 


i=l 


NA(t) 


N-l  N 

Z  £ 

i=l  j=i+l 


EB^t) 


EBj(t) 


NA( t)  NA( t) 


(14) 


when 


EB1(t) 

NA(t) 


«  1.0  , 


(15) 


when 


N 

EN(t)  =  £  EB , ( t )  =  EB( t)  (16) 

i=l 


Thus,  if  Q  =  1  for  all  components  and  E^(t)  is  very  small  (relative  to 
the  number  of  aircraft)  for  each  component,  then  EB(t)  is  an  approxima¬ 
tion  to  the  expected  number  of  non-mission-capable  aircraft  with  no  can 
nibalization  and  will  always  overstate  EN(t). 

U.  NMC  with  partial  cannibalization. 

Here  we  assume  that  some  components  are  relatively  easy  to  canni¬ 
balize,  and  that  some  are  so  difficult  to  remove  or  install  that  it  is 
not  desirable  to  cannibalize  them.  Let 


Ic  =  {  i  |  component  i  is  cannibalizable  ) 
In={  i|  component  i  is  not  cannibalizable  } 
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We  first  compute  the  probability  of  exactly  k  NMC  aircraft  due  to  items 
in  the  set  1^.  This  is  merely  the  PN^  (t)  derived  earlier  using  P(j), 
where  P(j)  is  computed  using  the  i  in  1^.  We  will  denote  this  probabil¬ 
ity  PNj(t).  The  probability  that  an  arbitrary  aircraft  is  nonopera- 
tional,  considering  only  shortages  of  cannibalizable  items,  is  then 


NA{t) 


j=l 


_j_ 

NA(  t) 


PNj(t) 


ENC(t) 
NA(  t) 


(17) 


The  probability  that  an  arbitrary  system  is  operational  after  shortages 
of  noncannibalizable  components  only  is 


II 

id 


(■ 


NA(  t )  J  ’ 


(18) 


assuming  that  Q^,  the  quantity  of  item  i  per  system,  is  one  for  each 

part,  i,  in  the  set  I  .  The  probability  that  an  arbitrary  aircraft  is 

n 

not  operational,  assuming  independence  of  demands  and  that  the  cannibal¬ 
ization  of  the  items  in  1^  takes  place  with  no  information  of  the  failed 
items  belonging  to  I  ,  is 


1  - 


ENC(t) 
NA(  t) 


II 

id 


(1 


EB^t) 

nauf  : 


1 


(19) 
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and  the  expected  number  of  nonoperational  aircraft,  given  partial  canni¬ 
balization  and  given  NA(t) ,  total  aircraft,  is 


ENP(t)  =  NA( t)  [  1  - 


1  - 


ENC(t) 
NA(  t) 


II 

iel 


1  - 


EB^t)' 

NA(t) 


(20) 


ENp(t)  is  an  overestimation  of  the  expected  NMC  if  information  describ¬ 
ing  which  noncannibalizable  items  had  failed  is  used  (since,  in  this 
case,  the  strategy  might  be  to  move  "holes"  due  to  cannibalizable  items 
to  those  aircraft  with  holes  due  to  noncannibalizable  items.  This  was 
not  considered  in  the  derivation). 

5.  Probability  of  meeting  aircraft  miss  ions . 

This  performance  measure  is  for  the  situation  in  which  it  is  desir¬ 
able  to  have  enough  operational  aircraft  to  perform  a  certain  demanded 
level  of  activity.  Let  B(t)  be  the  maximum  number  of  allowable  NMC  air¬ 
craft  that  still  permits  the  missions  to  be  met,  and  let  D(t)  be  the 
number  of  missions  at  time  t.  Let  r(t)  be  the  maximum  number  of  sorties 
per  unit  time  achievable  by  a  single  mission-capable  aircraft.  Note 
that 

B(t)  =  NA(  t)  -  ,  (21) 

where  [x]  is  the  ceiling  of  x,  that  is,  the  smallest  integer  so  that 
[x]  >  x,  and  D(t)  >  r(t).  NA(t)  is  assumed  to  be  an  integer.  The  probabil¬ 
ity  of  meeting  mission  demands  is  then  PP(t)  =  P(j)  (given  earlier) 
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where  j  —  B(t).  That  is,  the  probability  that  mission  demands  are  met 
is  simply  the  probability  that  the  number  of  nonoperational  aircraft  is 
less  than  or  equal  to  B(t). 

6.  The  expected  number  of  mission  demands  met ,  given  k  non- 
miss  ion  -capable  aircraft  is 


!D(t) 

r(t)(NA(t)  -  k) 


if  k  <  B(t) 

if  k  >  B(t) 


(22) 


Unconditioning  by  multiplying  by  the  probability  of  k  non-mission- 
capable  aircraft  (described  in  the  previous  paragraphs),  we  obtain 


NA(t) 

ES(t)  -  D( t)P( B( t) )  +  ^  r(t)(NA(t)  -  k)PNk(t)  ,  (23) 

k=B( t)+l 

where  P(B(t))  =  P(j),  j  =  B(t),  and  PNk(t)  were  described  in  paragraph  2 
above . 

7.  The  probability  distribution  of  the  number  of  mission  demands 


met . 

This  is  given  by 


PSk(t> 


K»(t)> 

II 

PN  (t) 
i  J 


if  k  -  D(t) 

if  k  -  r(t)(NA(t)  -  j); 

j  -  B(t)  +  1 . NA( t) 


(24) 


0 


otherwise  , 
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and  allows  the  variance  in  mission  demands  to  be  determined  from 

NA(t) 

VS( t)  =  D2 ( t)P(B( t) )  +  ^  (r(t) (NA(t)  -  k))2PNk(t)  -  ES2(t)  . 

k=B(t)+l  (25) 

Note  that  the  mission-related  measures  bring  into  play,  in  a  determinis¬ 
tic  sense,  additional  time-varying  parameters  such  as  the  number  of  air¬ 
craft,  NA(t),  the  demands  on  all  systems,  D(t),  and  the  portion  of 
demand  a  single  aircraft  can  satisfy,  r(t).  In  the  deployment  of  an 
aircraft  squadron,  NA(t)  might  represent  a  time-phased  deployment  of  the 
aircraft  as  well  as  attrition  due  to  losses  in  a  conflict.  The  demand 
D(t)  would  be  the  time-dependent  demand  for  sorties  as  an  engagement 
proceeds,  and  r(t)  would  be  time-varying  because  of  a  short-term  capa¬ 
bility  to  provide  surges  of  flying  activity. 
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V.  AN  INDENTURE  MODEL  FOR  TIME -DEPENDENT  PIPELINES 

Aircraft  components  are  classified  in  terms  of  major  assemblies, 
which  fit  directly  on  the  aircraft,  and  subassemblies,  which  fit  on 
other  assemblies.  In  the  Air  Force  the  major  assemblies  are  called  Line 
Replaceable  Units  or  LRUs,  and  the  subassemblies  are  called  Shop 
Replaceable  Units  (SRUs).  Further  breakdowns  of  the  assemblies  are  usu¬ 
ally  possible;  in  some  cases  this  breakdown,  or  indenturing,  goes 
through  five  or  more  levels.  Assemblies  and  their  subassemblies 
("Assembly"  and  "component"  will  be  used  interchangeably  in  this  sec¬ 
tion)  affect  the  parent  aircraft  in  different  ways.  A  shortage  of  an 
assembly  will  cause  a  hole  in  an  aircraft.  A  shortage  of  a  subassembly 
will  cause  a  hole  in  an  assembly,  which  may  or  may  not  cause  a  shortage 
of  the  assembly  and  a  consequent  aircraft  hole.  Thus,  in  terms  of 
effect  on  system  performance  and  in  terms  of  spare  parts,  subassemblies 
and  their  parent  assemblies  should  be  treated  differently. 

THE  CLASSICAL  STEADY-STATE  METHOD 

The  classical  method  of  treating  subassemblies  (see,  e.g.,  Muck- 
stadt,  1973)  is  based  on  a  number  of  steady-state  assumptions  that  do 
not  hold  true  for  the  time-dependent  case.  To  explain  the  dynamic  model 
it  will  be  useful  to  first  derive  the  steady-state  model  and  illustrate 
these  assumptions.  Let  d^  be  the  demands  per  unit  of  time  for  subcom¬ 
ponent  j  when  component  i  is  repaired.  Let  EB^  be  the  expected  back¬ 
orders  for  component  j  and  let  d_^  be  the  demands  for  repair  of  component  i 
per  unit  of  time.  The  average  repair  time  for  component  i  is  assumed  to 
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be  T\  .  The  first  step  in  the  classical  approach  is  to  determine  the 
average  waiting  time  for  subcomponent  j,  given  that  it  is  needed  in  the 
repair  of  i.  This  is  given  by  classical  queueing  theory  as 


W  .  =  — 

J  -5 


EB  , 


(1) 


EB^  represents  the  average  shortages  of  subassembly  j  that  must 
be  resupplied  before  new  demands  for  j  can  be  satisfied.  In  steady 
state,  the  rate  of  demand  fulfillment  is  equal  to  the  rate  of  demand,  so 
that  dj  represents  the  average  rate  at  which  backorders  for  j  are  satis¬ 
fied.  Then,  EB^/d^  gives  the  average  time  until  the  shortages  are 
cleared,  or  until  any  new  demand  can  be  satisfied. 

In  the  most  commonly  used  approach,  the  conditional  probability 
that  the  assembly  i  requires  subassembly  j  when  it  fails  is  approximated 
by  the  relative  demand  rates  for  j  and  i, 


d1 


ij 


(2) 


Under  the  assumption  that  subassembly  failures  are  independent,  the 
expected  time  that  assembly  i  will  wait  for  all  of  its  subassemblies  is 

EB  . 

_ J 

d  } 

subassemblies  j  ^ 


subassemblies  j 


(3) 


-42- 


This  time  is  then  added  to  the  assembly  repair  time  and  multiplied  by 
the  demand  rate  to  give  the  average  number  of  component  i  in  repair  plus 
awaiting  parts  (AWP) .  That  is, 

A  (in  steady  state)  =  •  (T  +  W^) 

■  diTi +  ebj  <4> 

subassemblies  j 

Steady-state  Sum  of 
pipeline  subassembly 

without  AWP  backorders 

This  result  requires  the  assumption  that  the  subcomponent  failure  sets 
are  nonintersecting  (which  really  means  that  no  more  than  one  subcom¬ 
ponent  can  fail  or  be  demanded  in  the  repair  of  each  assembly).  It  also 
requires  the  assumption  that  there  is  no  cannibalization  of  subassem¬ 
blies  (which  would  minimize  the  number  of  components  awaiting  parts). 
Finally,  it  assumes  that  a  subcomponent  does  not  appear  on  more  than  one 
type  of  parent  component. 


THE  MODEL  FOR  TIME-DEPENDENT  PIPELINES 

The  approach  used  for  dynamic  systems  is  similar  but  relaxes  some 
of  the  previous  assumptions.  Let  A^(t)  represent  the  number  of  com¬ 
ponents  of  type  i  awaiting  parts  at  time  t,  and  let  A^(t)  represent  the 
expected  value  of  A^(t).  The  mean  number  of  components  in  repair  and 
awaiting  parts  is  then[l] 


\(t)  =  Xt(t)  +  AA(t) 


(5) 


[1]  Let  Fi(t)  represent  the  probability  distribution  of  Aj(t)  and 
F^(t)  be  the  cumulative  Poisson  or  negative  binomial  distribution 


When  subassemblies  are  cannibalized,  EB^(t)  is  determined  in  a 
manner  similar  to  the  treatment  of  aircraft  cannibalization.  Let  P^(k 
be  the  probability  that  shortages  of  the  jth  subassembly  are  less  than 
or  equal  to  k.  Then 


Sj(t)+k 


PJ(k)  = 


Probability  of  l  failures  of  j 


from  the  Poisson  or  negative 
binomial  distributions 


described  earlier  for  these  components  in  repair.  When  Fj^(t)  and  F^(t) 
are  independent  and  both  Poisson,  or  when  they  are  both  negative  binomi¬ 
al  with  the  same  variance-to-mean  ratio,  the  distribution  of  the  com¬ 
ponents  awaiting  parts  and  in  repair  will  be  of  the  same  form  with  the 
mean  X^(t) .  Otherwise,  when  the. distributions  have  different  forms, 
convolutions  between  F^(t)  and  F^(t)  must  be  performed.  Although  it  is 
likely  that  F^(t)  differs  from  F^(t)  in  distributional  form,  our  experi¬ 
ence  has  shown  that  the  approximation  of  F^(t)  by  a  Poisson  or  negative 
binomial  distribution  does  not  overly  distort  the  results.  We  also  note 
that  convolution  to  obtain  the  resulting  distribution  is  usually  avoided 
in  the  steady-state  calculations  as  well. 
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Let  P^(k)  be  the  probability  that  the  number  of  nonoperational  assem¬ 
blies  of  type  i  (due  to  subassembly  shortages)  is  less  than  or  equal  to 
k.  Then,  since  this  probability  is  the  same  as  the  probability  that  all 
subassemblies  have  .  hortages  less  than  or  equal  to  k, 


Pt(k) 


n  pJ(k) 

subassemblies 
j  e  i 


The  expected  number  of  assemblies  of  type  i  that  are  unavailable  due  to 
subassembly  shortages  is  then  found  by  summing  the  tail  of  this  distri¬ 
bution: 


NA(t).Qi+Si(t) 


At(t)  - 


(1  -  PjCk))  , 


where  =  quantity  of  assembly  i  on  an  aircraft, 


NA(t)  =  number  of  aircraft  at  time  t,  and 


S  (t)  =  supply  level  for  component  i  at  time  t. 


The  distribution  of  parts  of  type  i  awaiting  subassemblies  is 


FA(t’k)  =  Pi(k)  “  Pi(k  "  1}  • 
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No  Cannibalization  with  Joint  Failures  Allowed 


»  \ 


An  alternative  form  of  EB^(t)  is  used  when  there  is  no  cannibaliza¬ 
tion  of  subcomponents.  The  equations  are  similar  to  those  used  for 
determining  the  number  of  aircraft  that  are  non-mission-capable  under  a 
policy  of  no  cannibalization.  In  this  case, 


Aj(t)  =  (NA( t)  .  Qi  +  Si(t)) 


1  - 


n  ii  -  — 

I  1  \  NA.Q  +S  (t)l 

qqpttjKI  lee  1  / 


subassemblies 


(11) 


When  there  is  more  than  one  subassembly  of  type  j  on  assembly  i,  this 
expected  value  takes  on  the  more  complicated  form  shown  in  Sec.  IV.  The 
probability  distribution  of  subassembly  shortages  in  the  no¬ 
cannibalization  case  is  not  known  although  its  moments  can  be  obtained. 
It  is  expedient,  although  not  precisely  correct,  to  use  the  geometric 
distribution  with  A^(t)  as  the  mean. 


Some  Cautions  and  Comments  on  the  Use  of  This  Model 

In  the  equations  used  for  determining  X_.(t),  one  must  be  careful  to 
avoid  overdrawing  the  influence  of  shortages  of  those  components.  For 
example,  when  no  repair  capability  exists  for  the  parent  assembly,  there 
can  be  no  await ing-parts  quantity  for  that  parent,  even  though  there  may 
have  been  subassembly  failures.  In  this  case  the  demand  for  the 
subassemblies  should  be  deferred  to  the  time  at  which  repair  occurs. 
Thus,  d  (t)  represents  the  discovered  failures  instead  of  true  failures 
of  the  subcomponents. 


Independence  between  subcomponent  and  parent  failure  distributions 
has  been  a  strong  assumption  in  much  of  this  derivation  (as  it  is  in  the 
steady-state  case),  even  though  most  parent  component  failures  are  due 
to  problems  with  the  subcomponents.  Clearly,  the  number  of  discovered 
subassembly  failures  cannot  exceed  the  number  of  failures  that  exist  on 
assemblies  that  have  entered  repair.  This  requires  conditioning  the 
distribution  of  the  subassembly  demands  on  the  number  of  parent  com¬ 
ponents  in  the  pipeline  and  determining  the  pipeline  and  AWP  quantities 
simultaneously.  This  means  that  the  foregoing  equations,  which  assumed 
independence,  will  sometimes  overstate  the  number  of  assemblies  in  the 
AWP  state  with  no  cannibalization  because  the  failed  subassemblies  are 
allowed  to  be  distributed  over  some  non-failed  assemblies.  It  also 
means  that  the  effects  of  subassembly  cannibalization  are  sometimes 
overstated,  since  the  cannibalization  is  allowed  to  take  place  across  a 
wider  range  of  assemblies  than  those  that  have  failed.  In  general,  how¬ 
ever,  the  independence  assumption  leads  to  reasonable  approximations 
since  the  rate  of  each  subassembly  failure  is  usually  considerably 
smaller  than  the  assembly  failure  rate  (due  to  each  assembly  usually 
having  large  numbers  of  subassemblies). 

Shared  Subassemblies 

For  simplicity,  the  development  in  this  section  was  based  on  an 
assumption  that  subassemblies  do  not  fit  on  more  than  one  type  of 
parent.  This  is  not  a  necessary  assumption  and,  when  violated,  it  can 
be  dealt  with  by  allocating  the  subassembly  shortage  distribution  back 
to  the  potential  parents  based  on  the  quantity  of  each  type  of  parent  in 
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the  aircraft,  Q^,,  the  quantity  of  the  subassembly  in  the  aircraft,  q., 

T 

and  the  quantity  of  the  subassembly  on  each  parent,  P .  Let  FA^(t,k) 
be  the  total  shortage  distribution  of  subassembly  j  for  the  aircraft. 


FA.  (t,k)  =  Probability  that  there  are  exactly  k  shortages 
of  component  j  at  time  t  . 


We  assume  here  that  the  failure  rate  of  j  is  not  affected  by  which 
assembly  it  is  on,  that  shortages  are  allocated  to  assemblies  based  on 
relative  demand,  that  assembly  i  has  P„  of  j-type  subcomponents,  that 
there  are  of  i-type  components  on  the  aircraft,  and  that  there  are  Qj 
of  j-type  components  on  the  aircraft.  Then  the  probability  that  there 
will  be  £  shortages  on  assembly  i,  when  there  are  k  (k  >  £)  total  shor¬ 
tages  of  subassembly  j,  is  given  by  the  hypergeometric  distribution, 


When  the  total  shortages  k  are  less  than  £,  there  is  of  course  0  proba¬ 
bility  of  1  shortages  on  assembly  i.  To  illustrate  this,  assume  =  £, 
P„  =  2,  and  =  6.  Then  2/6  of  the  subcomponents  j  lie  on  assembly  i. 
Given  4  shortages  of  j,  the  probabilities  of  1  and  2  shortages  on  i  are 


P(l,4) 


8 

15 


P(2,4) 
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The  probability  of  l  shortages  of  j  associated  with  i,  FAj(t,£), 
is  removed  from  its  dependence  on  the  number  of  shortages  of  j  by 
summing  across  the  possible  values  of  k.  That  is. 


In  case  only  the  expected  value  associated  with  FAj(t,£)  is  needed  (as 
is  the  case  in  the  no-cannibalization  assumption) ,  it  is  easy  to  show 


that 


"AS(t) 


E(FAj(t,*) 


EBj(° 


Thus,  the  total  expected  shortages  of  j  are  allocated  to  assembly  i 
based  on  the  number  of  subassemblies  on  i  relative  to  the  total  of  j  on 
the  aircraft.  The  subcomponent  shortage  allocation  described  in  the 
above  paragraph  is  performed  before  computing  the  effect  of  subcompo¬ 
nent  shortages  on  the  parents  under  the  cannibalization  and  no¬ 
cannibalization  assumptions.  Therefore,  there  is  an  implicit  assump¬ 
tion  of  no  cannibalization  of  subcomponents  across  unlike  parents. 


t 


-49. 

VI.  A  TIME -DEPENDENT  PIPELINE  MODEL  FOR  MULTIPLE 
ECHELONS  OF  REPAIR 


1 


* 


* 


Component  repair  is  performed  at  several  levels  within  the  Air 
Force.  Flight-line  personnel  remove  and  replace  major  assemblies. 
Squadrons  or  bases  repair  these  in  the  field  if  possible  or  ship  them  to 
larger,  centralized  facilities  serving  several  units.  This  section 
describes  a  mathematical  model  for  these  various  echelons  of  repair. 
Except  for  time  delays  for  transportation,  the  approach  is  very  similar 
to  that  used  to  model  indentured  components.  As  before,  the  commonly 
used  steady-state  approach  for  multiple  echelons  is  explained  first,  to 
form  a  basis  for  understanding  the  time-dependent  approach. 


THE  STEADY  STATE  MODEL 

The  decision  to  ship  a  component  to  another  echelon  is  based  on 

policy  (for  example,  perform  no  local  repair  on  this  component)  or  on 

random  chance  (inability  to  repair  certain  types  of  failures  locally). 

In  either  case  there  is  a  decision  to  ship  a  component  to  another 

echelon  of  repair.  At  the  time  a  component  is  shipped,  another  is 

ordered  (usually  from  the  same  echelon).  In  a  steady-state  model,  the 

average  rate  of  flow  to  the  next  echelon  from  a  unit  is  equaled  by  the 

average  rate  of  flow  back  to  the  unit.  Let  this  average  flow  for  unit  k 

F 

be  given  by  d^.  Assume  there  is  a  transportation  delay,  T^,  for  the 


retrograde  movement  (to  the  higher  echelon)  and  a  transportation  delay, 
F 

T^,  for  forward  shipments  (from  the  higher  echelon).  The  average 


repair  time  at  the  centralized  facility  is  given  by  x.  The  average 
number  of  components  in  repair  or  in  transportation  is  then  X  ,  given 

K. 


xk  =  5k(Tk  +  T) 


and  the  average  number  of  components  in  shipment  to  the  unit  is  X  , 

K- 

given  by 


Now,  if  there  is  only  one  unit  supported  by  this  higher  echelon,  the 
steady-state  model  is  quite  simple.  The  average  total  pipeline  for  a 
given  component  is  simply  the  base  repair  pipeline  for  the  component, 
XR,  plus  the  above  pipeline  averages: 


Xk  +  Xk  +  Xk  =  Ak  +  *k  (Tk  +  T  +  Tk) 


Given  a  total  stock  level  for  the  two-echelon  system,  the  various  per 
formance  measures  such  as  backorders  and  fill  rate  can  be  easily  com¬ 
puted  as  shown  earlier. 

Frequently,  the  total  system  stock  is  segregated  into  a  supply 
level  S  to  support  the  retrograde  and  repair  pipeline  of  the  higher 
echelon,  and  a  supply  level  to  support  the  base  repair  and  forward 
shipment  time  from  the  supporting  facility.  In  this  case  a  response 
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delay,  Ag,  is  computed  and  added  to  the  forward  transportation  time. 

The  response  delay  is  approximated  by  computing  the  backorders  in  the 

retrograde  and  repair  pipeline  when  the  supply  level  is  S  ,  These  are 
R  R 

denoted  EB(S  ,  Under  the  continuing  assumption  that  the  system  is 

in  steady  state,  these  backorders  or  shortages  (actually,  inability  to 
respond  to  orders)  are  cleared  at  the  rate  d^.  The  time  to  clear  all  of 
them,  and  therefore  the  time  delay  in  responding  to  a  new  order  (since 
the  backorders  represent  orders  already  placed),  is  given  by 


ERfSR  XRX 


(4) 


The  forward  pipeline  is  then  increased  by  this  time  delay  so  that 


xk  ’  5k  •  (Tk  +  V 


=  +  EB(SR,AR) 


(5) 


It  should  be  apparent  that  this  separation  of  supply  levels  does 
not  give  the  same  performance  as  that  computed  when  we  consider  one 
supply  level  for  the  entire  system.  That  is, 


Xk  = 


Xl  + 

k  k  k 


(b) 


with  supply  level  +  S  ,  is  not  the  same  as 

i  +  Vk  +  SBCS8,^)  (7) 

with  supply  level  S^.  In  the  first  case,  the  retrograde  and  repair 
stock  can  be  used  to  support  shortages  at  the  base  or  in  the  retro¬ 
grade  pipeline,  while  the  second  model  does  not  permit  the  retrograde 
supply  level  to  be  used  for  base  or  unit  support. 

If  only  one  base  were  supported,  it  would  always  make  sense  to 
place  all  stock  in  the  base  or  unit  supply  level  S^,  since  performance 
is  ultimately  measured  at  the  aircraft  and  not  in  the  retrograde  pipe¬ 
line.  When  this  is  the  case,  the  two  models  are  the  same,  since 


=  A®  +  dfcT^  +  EB(0,Ak) 


+ 


(8) 


When  multiple  units  or  locations  are  supported  by  higher-echelon 
repair,  it  is  necessary  to  consider  the  pipelines  from  each  location 
jointly  unless  there  is  ownership  of  specific  components.  Generally,  we 
assume  that  once  components  enter  the  higher  echelon  for  repair  they  can 
be  given  back  to  any  location.  It  is  difficult,  however,  to  determine 
performance  of  any  particular  unit  without  designating  a  supply  level 
for  it.  When  there  is  no  retrograde  supply  level,  the  unit  performance 
can  be  determined  without  considering  the  other  units  by  adding  the 
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retrograde  and  repair  pipelines  of  each  individual  unit  and  computing 
measures  using  the  unit  supply  levels.  The  same  performance  for  the 
system  as  a  whole  can  usually  be  achieved  at  lower  supply  cost  by  pro¬ 
viding  a  single  retrograde  plus  central  repair  supply  level  for  all 
retrograde  pipelines.  In  this  case  we  compute  a  delay,  as  shown  ear¬ 
lier,  using 


X 


R 


£  * 

k  bases 


E 

k  bases 


5k(T£  + 


(9) 


and  then  computing 


A 


S 


EB(SR,XR) 


_  y 


k  bases 


(10) 


which  is  the  delay  in  responding  to  any  unit's  order.  The  unit  perfor¬ 
mance  is  then  based  on 


k 


VTk + 


V 


Ak  +  dkTk  +  eb<sR>xR) 


(id 


k  bases 
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An  alternative  explanation  of  this  model  is  that  the  retrograde  repair 
R  R 

shortages,  EB(S  ,  X  ),  are  allocated  on  the  basis  of  the  relative  aver¬ 
age  demands  placed  by  the  units. 


TIME -DEPENDENT  MODEL 

The  dynamic  multiechelon  model  requires  introduction  of  the  follow¬ 
ing  variables: 


X^(t)  Average  number  of  components  from  base  k  in 

transportation  to  or  in  repair  at  the  next  higher 
echelon  at  time  t. 


x“(t) 


Average  number  of  components  in  base  repair  at 
time  t. 


x[(t) 


Average  number  of  components  in  transportation 
to  base  k  at  time  t. 


sR(t) 

sk(t) 


Supply  level  for  the  retrograde  pipelines  at 
time  t. 

Supply  level  at  base  k  at  time  t. 


The  average  shortage  at  time  t  in  the  retrograde  transportation  pipeline 
is  given  by  the  expected  backorders , 


EBR(SR(t),AR(t))  . 


Now,  if  location  k  is  the  only  one  supported,  then  the  shortages  seen  at 
that  location  are  given  by  the  excess  of  local  repairs  and  unfulfilled 


I 


orders  (placed  with  the  higher  echelon)  over  the  local  supply  level, 
S^(t).  The  average  number  of  unfulfilled  orders  at  time  t  is  given  by 


Xjj(t)  =  Xj(t)  +  EBR(SR(t),  XR(t))  . 


Depending  on  the  repair  function  at  the  higher  echelon,  the  average 
number  of  components  in  the  forward  pipelines  can  be  difficult  to  deter¬ 
mine,  because  the  number  must  be  conditional  on  the  completion  of  repair 
and  retrograde  transportation.  A  useful  alternative  representation  of 
the  unfulfilled  orders  is  ,  ven  by  the  demands  placed  during  one  forward 

transportation  time  (since  all  of  these  must  be  unfulfilled)  plus  the 

F 

shortages  that  existed  at  time  t  -  T  in  the  retrograde  plus  repair  por¬ 
tion  of  the  system.  The  average  unfulfilled  orders  are  then 


.u.  .  /*  ,R,  ,,  .  „F.  ,R.  F 

-  /  F  dk(s)ds  +  EB  (S  (t  -  Tk),Xk(t  -  Tfc 

Jt-T, 


and  the  average  pipeline  at  location  k  is 


*k(t) 


dR(s)ds  +  EBR(SR(t 


V’xk(t  -  Tk}) 


+ 


X|j(t)  +  X®(  t)  . 


The  development  for  multiple  locations  or  squadrons  served  by  one 


higher  echelon  repair  facility  is  similar,  except  we  assume  that  the 
supply  level,  SR(t),  applies  to  all  retrograde  pipelines  and  that  the 
shortages  must  be  allocated  back  to  the  units.  Thus, 


At)  =  2  x*(o 

k  bases 


(16) 


and  the  average  shortages  in  the  retrograde  portion  of  the  system  are 


EBR(SR(t),AR(t))  . 


(17) 


These  shortages  should  be  allocated  to  the  various  locations  or  squad¬ 
rons  based  on  a  criterion  representative  of  shortage  allocations  in 
real  situations.  One  criterion  suggested  by  steady-state  methods  is  the 
relative  number  of  demands  placed  at  time  t.  That  is,  the  shortages 
allocated  to  location  k  at  time  t  are  given  by 


<(t) 


<(t) 


EBR(SR(t),  \R(t)) 


k  bases 


(18) 


This  represents  a  rather  volatile  reallocation  of  shortages  when  the 
demands  at  the  various  locations  go  through  sudden  relative  changes.  An 
alternative  allocation  is  based  on  the  time-averaged  demand  at  the  loca- 


tion.  Let  D^(t,x)  be  the  average  demands  in  the  interval  (t  -  t,  t) . 
That  is. 


D^(t.t) 


d*(t)dt 


(19) 


Then,  we  can  use  the  allocation 


EBk(t) 


D*(t,t) 


DE(t,T) 


♦  EBR(SR(t),  XR(t)) 


k  locations 


(20) 


The  larger  the  value  of  t,  the  more  the  shortage  allocation  will  be 
influenced  by  (or  averaged  across)  past  orders.  The  appropriate  value 
of  T  Is  still  open  to  empirical  testing,  but  it  seems  reasonable  to  use 
a  value  of  t  near  the  average  time  it  would  take  the  central  facility  to 
reallocate  shortages.  If  the  only  means  of  reallocating  shortages  by 
the  central  facility  is  to  not  fulfill  orders  selectively,  then  the  rate 
at  which  they  can  be  allocated  to  any  location  is  limited  by  the  rate  of 
orders  placed  by  that  location.  Note  that  when  the  reallocation  process 
is  limited  in  this  way,  the  shortages  allocated  to  any  location  are  lim¬ 
ited  by  the  number  of  orders  placed  by  that  location.  (That  is,  the 
unfulfilled  orders  are  limited  by  the  orders  placed.)  The  allocation  of 
backorders  using  the  above  equation  and  x  =  t  will  keep  the  average 
number  of  backorders  allocated  below  the  average  total  demands  from  0  to 


t,  since 


8°kV> 


oJ(t.O) 
k  locations 


B0R(SR(t),  XR( t) ) 


D*(t,0) 


X)  Dk(t’0) 

k  locations 


(21) 


X)  D*(t,0)  =  Dk(t>°> 

k  locations 


In  our  use  of  the  model,  we  have  approximated  the  distribution  of  allo¬ 
cated  shortages  with  the  Poisson  or  negative  binomial  distributions  with 
£ 

mean  value  BO  (t)  .  This  approach  is  consistent  with  the  steady-state 

K. 

approaches  commonly  in  use.  (See,  for  example,  Sherbrooke,  1968(a).) 

The  above  equations  describe  with  sufficient  generality  the  multi¬ 
ple  echelon  time -dependent  model.  Interruption  of  transportation  and 
subassembly  interaction  with  higher  echelons  are  additional  aspects  of 
this  model. 

INTERRUPTION  OF  TRANSPORTATION 

This  section  describes  a  model  that  has  fixed  or  deterministic 
transportation  times.  Nondeterministic  transportation  times  are  possi¬ 
ble  within  this  model,  but  the  description  of  unfulfilled  orders  is  con¬ 
siderably  more  difficult.  One  variation  that  is  easily  accommodated  is 
to  allow  interruptions  of  transportation  for  fixed  periods  of  time.  For 
interruption  of  the  retrograde  transportation,  the  pipeline  calculation 
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of  A^Ct)  is  modified  accordingly.  (See  the  earlier  section  regarding 

pipeline  calculations.)  The  cutoff  of  forward  transportation  affects 

u  F 

the  equation  for  X^(t).  If  the  cutoff  ends  before  time  t  -  T  ,  or 

F 

starts  after  t,  there  is  no  modification.  If  it  starts  before  t  - 
(at  S)  and  ends  after  t,  we  use 


x“(t)  -y  d£(s)ds  +  B0R(SR(S),  A*(S))  •  (22) 

F  F 

If  it  starts  prior  to  t  -  and  ends  in  the  interval  [t  -  T^,t],  we  use 

F 

this  also.  If  the  cutoff  starts  in  the  interval  [t  -  T^,t]  and  ends 
after  t,  then  the  original  no  cutoff  equation  applies. 

INDENTURED  COMPONENTS  AND  MULTIPLE  ECHELONS 

The  previous  section  discussed  the  approach  for  indentured  com¬ 
ponents.  In  that  approach  it  was  assumed  that  the  demands  for  subassem¬ 
blies  occurred  at  the  same  time  that  demands  for  repairs  of  the  parent 
assemblies  occurred.  This  is  a  reasonable  approximation  for  local 
repair,  but  is  not  correct  for  the  repair  of  subassemblies  at  a  facility 
remote  from  the  aircraft.  In  this  case  the  demand  for  repair  of  the 
major  assembly  occurs  immediately,  but  the  demand  for  a  subcomponent  is 
not  known  until  the  parent  component  is  received  and  tested  at  the 
remote  facility.  The  approach  taken  in  this  case  (for  higher  echelons) 
is  to  delay  the  demand  for  subcomponents  by  the  transportation  time  (or 
cutoff  time  plus  transportation  time  if  the  retrograde  transportation  is 
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cut  off).  Thus,  the  average  demand  from  location  k  for  a  subcomponent 
at  the  remote  facility  is  given  by 


if  no  cutoff  or  if 
cutoff  ends  before 
t  -  t^  or  starts 
after  t 


R 


if  t  falls  in  a  cut¬ 
off  period  or  cutoff 
ends  less  than 
days  prior  to  t 


if  t  -  Tfc  =  B  where 
a  and  B  indicate 
the  start  and  end  of 
retrograde  cutoff, 
respectively. 


With  this  definition  of  demand,  the  indenture  model  of  the  previous  sec¬ 
tion  operates  as  described  earlier. 


(23) 
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VII.  TIME -DEPENDENT  OPTIMAL  DETERMINATION  OF  SPARE  PARTS 
TO  MEET  AN  OPERATIONAL  OBJECTIVE 


Although  built  primarily  as  a  readiness  assessment  tool,  the  Dyna- 
METRIC  model,  because  of  its  analytic  framework,  permits  the  determina¬ 
tion  of  spare  parts  required  to  satisfy  a  given  level  of  aircraft  avail¬ 
ability.  In  fact,  this  enhances  the  readiness  assessment  by  giving  an 
answer  to  the  question  of  how  many  additional  spare  parts  are  needed  to 
provide  a  given  level  of  readiness  when  goals  are  not  achieved.  (Of 
course,  spare-parts  provisioning  is  only  one  of  several  ways  to  achieve 
a  given  level  of  capability.)  The  determination  of  supply  levels  within 
Dyna-METRIC  is  separated  into  four  phases,  dealing  with:  spare  parts  to 
overcome  queueing  in  test  facilities,  spare  parts  for  higher  echelons, 
spare  parts  for  subassemblies,  and  spare  parts  of  major  assemblies  at 
the  squadron.  Each  of  the  latter  three  phases  will  be  discussed 
separately  in  this  section.  The  spare-parts  determination  for  queueing 
is  described  in  Sec.  VIII.  First,  however,  the  general  approach  for  a 
time-dependent  stockage  policy  will  be  discussed. 

The  fact  that  pipelines  have  time-dependent  probability  distribu¬ 
tions  means  that  the  optimal  mix  of  spare  components  at  one  point  in 
time  may  not  be  the  optimal  mix  at  another.  For  example,  when  there  is 
no  repair  capability  for  some  period  of  time  for  a  subset  of  the  com¬ 
ponents,  the  correct  mix  of  spares  will  require  substantial  numbers  of 
these  components;  fewer  will  be  required  if  repair  capability  is  pro¬ 
vided  at  some  later  point  in  time.  Clearly,  there  is  a  mix  of  com¬ 
ponents  that  will  provide  the  desired  capability  at  each  point  in  time 


at  the  lowest  cost.  This  model  does  not  determine  this  cross-time 


optimal  value  for  several  reasons.  First,  the  time  horizon  is  fre¬ 
quently  a  policy  variable,  and  the  cost  of  the  spare  parts  is  desired 
for  marginal  increases  in  the  horizon.  That  is,  the  decisionmakers  gen¬ 
erally  want  to  know  the  incremental  costs  and  spare  parts  to  achieve  a 
given  capability  for  a  slightly  lengthened  time  horizon.  Complete 
changes  of  the  spare-parts  mix  as  the  time  horizon  increases  may  not  be 
desirable,  particularly  when  supply  levels  already  exist.  Practically, 
furthermore,  the  methodology  for  cross-time  optimization  has  not  been 
developed.  The  approach  taken  is  to  compute,  for  each  time  of  interest, 
the  marginal  increase  in  spare  parts  to  achieve  a  given  capability  over 
those  already  input  or  determined  for  a  previous  time .  Thus ,  embedded 
in  each  subsequent  time  is  a  supply  level  that  will  achieve  the  desired 
capability  at  all  previous  times. 

SUPPLY  LEVELS  FOR  SQUADRONS  OR  OPERATING  UNITS 

Aircraft  mission  capability  at  the  squadron  is  the  primary  interest 
and  is  the  constraint  in  determining  the  supply  level.  This  is 
expressed  as  a  probability  of  not  exceeding  a  given  number  of  non- 
mission-capable  aircraft  at  each  time  point  of  interest.  That  is,  the 
decisionmaker  selects  the  confidence  of  achieving  a  given  mission- 
capable  rate  for  the  squadrons.  In  determining  the  supply  level,  the 
model  attempts  to  provide  enough  spare  parts  to  give  the  desired  confi¬ 
dence  at  lowest  cost  at  each  point  in  time  of  interest.  Thus,  the 
objective  function  is  the  total  cost  of  spare  parts  at  the  squadron. 

Let  represent  the  spare-parts  level  for  component  i,  c^  the  unit  cost 
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of  component  i  (per-unit  costs  are  assumed  to  be  independent  of  quan¬ 
tity)  ,  a  the  desired  confidence  level,  the  non-mission-capable  rate 
not  to  be  exceeded,  and  P(K^,S)  the  probability  that  the  non-mission- 
capable  rate  is  less  than  given  a  stock  level  S  (the  vector  of  S^'s). 
Then  the  problem  solved  for  each  squadron  or  operating  location  at  each 
time  of  interest  is 


minimize 


i  components 


CiSi 


subject  to  P(K^,S)  >  a 


S,  >  S.  for  each  i, 

i  —  io 


Si  integer 


(1) 


where  is  the  input  stock  level  or  previous  time  optimization  stock 
level  for  component  i.  The  probability  function  P(KN»S)  is  the  form 
given  in  Sec.  IV  for  full  cannibalization.  That  is, 


P(Kn»S) 


n  pi(q1  • 

i  components 


V 


(2) 


where 

W*K 

P1(Q1  •  K^)  =  ^  P^k)  (3) 

k-0 

and  P^(k)  is  the  probability  of  exactly  k  failures  of  component  i. 
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The  reader  may  observe  that  a  necessary  condition  for  the  perfor¬ 
mance  constraint  to  be  met  is  to  have 


Pi(Q1.KN)  >  a 


(4) 


for  each  i  (since  otherwise  the  product  of  these  functions,  which  have 
values  less  than  or  equal  to  1,  would  not  be  greater  than  a). 

Thus,  the  first  step  in  the  requirements  process  is  to  find  a  stock 
level  for  each  component  that  satisfies  this  necessary  condition 
independently.  Once  this  step  has  been  completed,  the  programming  con¬ 
straint  will  not  usually  be  satisfied,  and  this  is  the  point  at  which 
optimization  comes  into  play.  Marginal  analysis  is  used  to  determine 
the  "best"  mix  of  additional  components  to  achieve  the  desired  goal. 

This  process  proceeds  by  investing  in  one  additional  component  at  a  time 
which  is  selected  by  finding  the  component  that  gives  the  largest 
increase  in  the  logarithm  of  the  confidence  level  at  the  lowest  cost.[l] 
That  is,  we  determine 


A^n  PO^.S) 


(5) 


[1]  The  logarithm  of  P ( Kjq , S )  allows  the  constraint  to  be  converted 
from  a  single  product  to  a  sum  of  terms.  This  separation  of  terms  then 
allows  marginal  analysis  to  find  an  optimal  or  nearly  optimal  solution. 
See  App.  A. 


r 


\ 
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where 


Ai  *n  P(Kn*S)  =  *n  I^<km.s1)/p(kn,s)] 


(6) 


and 


<S1*V 


,S,  +  1,  S 


i+1  ’ 


...) 


(7) 


The  component  for  which  supply  is  increased  one  unit  is  the  one  whose 
index  solves 


in  P(VS) 

max  - 

i  Ci 


(8) 


This  process  continues  until  the  given  confidence  level  is  achieved.  At 
this  point  the  resulting  supply  vector  S*  will  be  an  "undominated"  or 
"efficient"  solution  of  the  base  stockage  problem  (see  Fox,  1966). [2] 

For  most  purposes  an  undominated  solution  can  be  considered  to  be 
optimal  or  nearly  optimal. 


[2]  To  reach  the  undominated  solution,  the  values  of  a  and  KN  must 
be  set  to  "reasonable"  levels.  See  App.  A  for  more  discussion  of  this. 
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SUPPLY  LEVEL  FOR  INDENTURED  COMPONENTS 

Supply  levels  for  subcomponents  are  determined  by  a  heuristic  that 
attempts  to  weight  the  subcomponent  need  based  on  the  relative  value  of 
shortages  of  parent  assemblies.  At  the  same  time,  it  attempts  to 
achieve  a  minimum-cost  solution  among  the  subcomponents  and  parent 
assemblies.  Muckstadt  (1980),  in  the  mod-METRIC  model,  used  a  search  of 
the  budget  space  to  determine  an  optimal  solution  for  the  case  in  which 
system  performance  in  terms  of  supply  backorders  was  to  be  optimized. 
Kotkin  (1978)  described  efficient  cost  heuristics  for  the  dual  problem, 
which  minimized  cost  subject  to  a  backorder  constraint.  This  approach 
uses  similar  heuristics  for  determining  quantities  of  subassemblies  that 
minimize  cost  subject  to  the  operationally  ready  rate  measure  described 
in  the  previous  subsection. 

Recall  that  the  problem  of  interest  is  to  minimize 

Z'i  St 


subject  to: 


np] 

i 


WiVV  > a 


(9) 


S.  >  S.  for  all  i  . 
i  -  io 
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The  indenture  problem  adds  the  variable  Sj ,  the  supply  level  of  the 
subassembly  j.  In  this  case  the  problem  becomes: 


minimize 


?c‘s‘  +^c‘  SJ 


subject  to  n  Pi(Qi.KN,S1,SSi)  >  a 


>  0  for  all  i 


S®  >  0  for  all  j  , 
s  i 

where  S  represents  the  supply  of  all  subcomponents  that  have  component 

i  as  a  parent.  Note  that  the  variable  name  has  been 

s  i 

expanded  to  show  the  dependence  on  S  by  writing  it  as 


P1(QIKs.s1,SSl) 


With  the  logarithmic  transformation  this  becomes 


minimize 


£ctst  +^'jsSj 


X>n  P1(Q1KN,Si,SSl)  Hna 


St  >  0  for  all  i 


Sj  >  0  for  all  j 


subject  to 


(11) 
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An  appropriate  Lagrange  multiplier,  Y,  allows  the  constraint  to  be  moved 
to  the  objective. [3]  This  converts  problem  (11)  to  the  following  form: 

minimize  £c  S  +£C!S4  -  Y  (]£  *n  P1(Q 
i  j  2  J  \i 

subject  to  >  0  for  all  i  (12) 


si. 


iKN’Si’S  )-£n  “ 


Sj  >  0  for  all  j  . 


With  an  appropriate  value  of  Y  (we  will  show  how  to  find  an  approx¬ 
imate  Y  shortly),  this  problem  is  equivalent  to  (11).  When  Y  and  are 

g 

known,  it  is  relatively  easy  to  determine  for  each  j.  The  overall 
heuristic  is  to: 


s  so 

1.  Start  with  a  given  value  of  S.,  say  S.  . 

3  3 


2.  Determine  under  the  procedure  given  earlier  in  this 
section. 

3.  Determine  an  appropriate  value  of  Y. 


4.  Determine  S^.  given  Y. 


5.  Redetermine  S.  given  the  new  values  of  S.. 

i  J 


[3]  Fox  (1978)  and  Fox  and  Landi  (1969)  discuss  the  existence  of 
this  multiplier.  It  can  be  thought  of  as  a  cost  or  penalty  for  not 
satisfying  the  performance  constraint. 
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Once  step  2  has  been  completed,  the  value  of  Y  can  be  approximately 

g 

determined  by  assuming  that  S_.  is  constant  and  applying  an  optimality 
condition  to  problem  (12).  [4]  That  is,  if  S  =  (S.^  Sfc)  solves 

(12),  then 


Ci(Si  +  -  ci(Si)  "  Y  t£n  pi(QiKN,Si  +  1,SS1)  "  £n  pi^iKN*Si»sSi)]  ^  ° 


for  all  i 


(13) 


« 


and 


Ci(si  ~  ~  W  ■ Y  t£n  pi(QiVsi  -  i.ssi) 


(14) 


-  in  Pi(QiKN,Si>SSl)]  >  0 
for  all  i 


« 

Thus , 


Y  < 


in  [P1(Q1Kn.S1  +  l,Ssi)/Pi(QiKN,S1,5si)] 


(15) 


[4]  This  condition  merely  states  that  the  optimal  value  of  the  vec¬ 
tor  S  must  give  a  smaller  value  of  the  objective  than  does  any  other 
value . 


* 
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and 


Y  > 


Jin  [Pi(QiKN,Si,SSl)/Pi(Q1KN,Si  -  l,SSi>] 


(16) 


which  provide  bounds  on  the  value  of  Y.  The  algorithm  uses  the  upper 
bound,  which  will  err  on  the  side  of  too  many  subassemblies  (which  typi¬ 
cally  cost  much  less  than  the  parent  assemblies).  That  is,  it  uses 


min 


i  in  [Pi(Q1KN,S1  +  l,S8l)/P1(Q1KN,S1,  SSl)] 


(17) 


This  is  particularly  convenient  since  it  is  the  reciprocal  of  the  ratio 
used  in  the  solving  (1)  at  step  2,  and  the  index  solving  this 
corresponds  to  the  last  component  of  supply  increased  in  that  process. 
(Thus,  the  approximate  Y  is  immediately  known  because  it  is  the  marginal 
improvement  ratio  for  the  last  stock  increase.) 

Given  this  value  of  Y  the  subassembly  problem  is  next  solved  in 
step  4  of  the  heuristic.  That  problem  is,  given  primary  assembly  stock, 


minimize 


j  J  J 


-  Y  in  P1 


<QiVsi»sSl) 


-  in  cx)j 


(18) 


S®  >  0 


for  all  j  . 
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Appendix  B  shows  that  this  problem  is  solved  approximately  by  find¬ 
ing  a  level  of  stock  for  each  subcomponent  j  that  satisfies 


“j  <v  * 1 


where  RRj(S^)  is  the  .ready  rate  for  subcomponent  j  defined  in  Sec.  Ill 
and  Bj  is  a  constant  defined  in  App.  B. 


SUPPLY  LEVELS  FOR  HIGHER  ECHELONS 

The  description  of  multiple  echelons  in  Sec.  VI  indicated  that,  for 
one  base  or  squadron  served  by  one  higher  echelon,  it  was  more  efficient 
to  determine  one  system  supply  level  rather  than  provide  supply  levels 
fof  the  squadron  and  higher  echelon  separately.  This  subsection 
addresses  the  problem  of  providing  spares  to  a  higher  echelon  facility 
serving  several  squadrons  at  different  locations.  In  this  case  it  is 
usually  efficient  to  provide  a  system  supply  level  for  the  retrograde 
plus  repair  pipeline  of  the  higher-echelon  facility,  and  separate  local 
supply  levels  for  each  of  the  squadrons  to  cover  the  order  and  shipping 
time  of  ordered  components  and  for  those  repairs  done  at  the  squadron. 
For  the  centralized  supply,  we  employ  a  heuristic  similar  to  that  used 
for  indentured  components : 


Start  with  a  given  level  of  supply  at  the  higher  echelon,  say 


Determine  S.  at  the  squadrons  as  discussed  earlier  (after  com- 

CO 

puting  the  e r feet  of  central  system  shortages  using  ). 


3.  Determine  an  appropriate  set  of  Lagrange  multipliers  given  the 


solution  in  step  2. 

c 

4.  Determine  the  supply  level,  S^,  at  the  higher  echelon  given  the 
multipliers  and  the  from  step  2. 

5.  Redetermine  at  each  squadron  given  the  new  value  of  central 
system  supply. 


The  overall  problem  to  be  solved  is [5] 

minimize  J^c.S 

i  i  1 


subject  to  11  pJ(Q.K  S^.sJ)  >  a  for  all  k  locations  (20) 

served  [6] 


S^,S^  >  0  for  all  k  and  1. 

Using  Lagrange  multipliers  and  the  logarithmic  transformations,  this 
problem  becomes 

minimize  ]£  c  S*  +  S<T 

k  locations  i  i 

2  Y,(l>n  pJ(Q,K  isfj-Jna) 

k  locations  \  i  1  1  N  1  1  ' 

k  c  (21) 

S^,S^  >0  for  all  k  and  i. 


[5]  The  new  notation  indicates  the  dependence  of  the  constraint  on 
the  central  system  supply  rather  than  the  subcomponents,  and  also  shows 
the  dependence  on  a  location  dependent  supply  level,  S^. 

[6]  The  probability  expression  in  the  constraint  is  the  same  as 
that  in  problem  (1)  and  problem  (2)  of  this  section  and  described  ear¬ 
lier  in  Sec.  IV. 
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C  CO 

With  the  central  system  supplies,  S^,  held  constant  (at  ),  this 
decomposes  into  subproblems  involving  only  the  component  i  at  location  k 
and  the  squadron  or  base  supply: 


min 


CiSi  +  Yk(£n  PJ(QiKN>Si’Si°>  ‘  Zn  “> 


si  >  0  * 


(22) 


The  value  of  is  available  after  this  problem  is  solved  through  margi¬ 
nal  analysis,  as  shown  in  the  last  subsection.  (In  fact,  it  is  the  same 
Lagrange  multiplier  that  is  used  in  the  indenture  problem  for  the  squa¬ 
dron  or  base.)  At  step  4  the  problem  solved  is 


min 


:‘s‘ '  ? 


£n  pi(QiVsi»si) 


(23) 


S 


c 

i 


>  0  . 


If 

With  the  base-level  supply,  S^,  held  constant  at  the  value  given  in 

step  2 ,  this  problem  is  solved  by  increasing  the  central-system  spares 
c 

level,  S^,  until  the  objective  function  no  longer  decreases. 

The  reader  may  observe  that  since  the  model  of  the  higher  echelon, 
described  in  Sec.  VI,  actually  has  that  echelon  affecting  the  squadron's 
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r 


V 


performance  after  a  time  delay  equal  to  the  forward  transportation  time, 

the  point  in  time  for  determination  of  a  supply  level  is  ambiguous  for 

several  locations  with  different  transportation  times.  The  function 
Ic  Ic  c 

Pi  ^i’^N’^i’^i^ ’  ^°r  eac^  l°cati°n  accommodates  these  separate  times,  and 

c 

it  is  necessary  only  to  assume  that  remains  constant  during  the  dif¬ 
ferent  delay  periods. 

The  calculation  of  subassembly  supply  levels  at  the  higher  echelon 

c 

interacts  with  the  calculation  of  S^.  Given  the  subproblem  (23),  we  can 
break  this  into  one  involving  identical  components: 

min  c^J  -  I  Yk  in  pj(Q.K  ,S^,SJ,S®C) 
k 


subject  to 


*»r  ^ 


(24) 


where 


,sc 


(s®c,s8C, 


s8c) 

bk  '  » 


a  vector  of  subcomponents  associated  with  component  i  at  the  higher 
echelon.  Actually,  since  the  subcomponents  may  affect  more  than  one 
major  assembly,  we  need  to  solve 


mlr,  E  ctsj  -  E  E  t„  P^Q^.sJ.sJ.s”) 


sj.s82  >  0 


(25) 


for  all  i. 
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This  problem  is  made  separable  in  terms  of  the  subcomponents  by  holding 

SC  c 

and  fixed,  and  differentiating  with  respect  to  the  subcomponent 
shortage,  (Ssc) ,  as  shown  earlier.  We  then  obtain,  as  illustrated  in 
App.  B,  subproblems  with  the  form 


,  s_s  V*  „k  _  sc. 

min  Vl' 


(26) 


S 


where 


aP±(QiKN>S±»S±»S±C) 

3  B  (SSC) 

3  3  (27) 

which 

have 

subassembly  i 


components 


pik«A-sik^'S“> 


The  subproblem  can  be  easily  solved  by  increasing  the  value  of  until 
the  objective  function  no  longer  decreases.  The  complete  algorithm  is: 

1.  Start  with  a  given  level  of  supply  at  the  higher  echelon,  say 

CO  SCO 

and  S  ^  for  all  components  and  subcomponents. 

2.  Determine  at  each  of  the  k  squadrons  (after  computing  the 

CO  S  CO 

effect  of  central  system  shortage  using  and  S  ). 

3.  Determine  an  appropriate  set  of  Lagrange  multipliers  given  the 
solution  in  step  2. 
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4.  Given  the  value  of  S,  from  step  2,  and  the  initial  value  of  S. 

k  J 

k 

for  all  subcomponents,  determine  the  partial  derivative  6^  with 
respect  to  each  subcomponent  for  each  location. 


SC 

5.  Determine  the  optimal  value  of  for  the  central  system. 

c 

6.  Given  this  value,  determine  the  optimal  value  of  for  the 

central,  system  using  Yk  and  subproblem  23. 

esc  k 

7 .  Given  S .  and  S .  ,  redetermine  S . . 

i  J  i 


« 
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VIII.  LIMITED  SERVICE  CAPACITY 


A  very  basic  assumption  in  Hillestad  and  Carrillo  (1980)  and  in  the 
previous  sections  of  this  report  is  that  sufficient  slack  service  capac¬ 
ity  exists  to  avoid  queueing  in  the  repair  of  components.  Pokress 
(1977)  has  shown  that  the  infinite-server  assumption  is  valid  as  long  as 
average  demands  remain  less  than  80  percent  of  service  capacity.  But 
surges  of  aircraft  flying  may  cause  increased  failures  such  that  the 
component  repair  resources  (manpower,  facilities,  or  test  equipment) 
become  overloaded.  When  this  overload  goes  on  indefinitely,  the  repar¬ 
able  component  queues  can  be  very  large,  but  when  there  is  slack  capac¬ 
ity  after  a  surge,  the  queue  can  be  worked  down.  We  are  interested  in 
the  time-dependent  mission  capability  of  aircraft,  that  is,  how  fast  the 
mission  capability  degrades  and  recovers  from  such  variations  in  demand. 
The  steady-state  results  of  queueing  theory  for  limited-server  queues  do 
not  provide  information  about  this  transient  behavior,  and  the  few 
analytical  solutions  to  limited-server  queues  in  the  transient  case  are 
too  restrictive  to  capture  the  complexities  of  multiple-use  test  equip¬ 
ment,  priority  scheduling,  etc. 

Dyna-METRIC  approximates  the  queueing  probability  distribution  with 
a  simple  simulation.  The  basic  approach  has  been  to  decompose  the 
queueing  and  service  portion  of  the  pipeline.  That  is,  the  total  aver¬ 
age  pipeline  of  a  given  component  is  assumed  to  be  composed  of 


A(t)  =  A(t)  +  B(t)  , 


the  simulation  to  determine  queue  sizes,  it  is  necessary  to  kno  the 
probability  that  a  given  test  stand  is  available  at  a  random  point  in 
time.  The  following  assumptions,  based  on  empirical  studies  of  certain 
types  of  test  equipment  by  Gebman,  Shulman,  and  Moore  (forthcoming) 
allow  us  to  use  a  dynamic  pipeline  model  similar  to  those  in  Sec.  II  to 
determine  this  probability.  These  are: 

1.  The  rate  at  which  test-stand  part-shortages  are  generated  at 

T 

any  point  in  time  is  known  to  be  d  (t).  (Note  that  only  the 
overall  rate  of  shortage  is  represented,  and  that  particular 
components  of  the  test  stands  are  not  distinguished  in  this 
model . ) 

2.  The  resupply  for  test-stand  shortages  is  given  by  the  probabil- 

-T 

ity  function  F  (t,s),  the  probability  that  a  shortage  at  time 
s  remains  at  time  t. 

3.  Shortages  of  test-stand  components  are  consolidated  on  the 
smallest  number  of  stands  possible. 

4.  A  shortage  of  a  test -stand  component  has  an  associated  proba¬ 
bility,  p^,  that  an  aircraft  component  i  can  still  be  tested 
despite  that  shortage. 

5.  Shortages  occur  on  no  more  than  one  of  any  test-stand  component 
type  at  a  time,  so  that  the  number  of  test  stands  not  opera¬ 
tional  after  cannibalization  as  never  more  than  one  stand. 

Using  these  assumptions,  we  can  determine  the  average  number  of  test- 
stand  shortages  with 


ffl 


The  probability  of  having  exactly  k  shortages  is  then,  under  a  Poisson 
rate  of  shortage  generation, 


PT(t,k) 


xT(t)k.-»T<t) 

k! 


(2) 


The  probability  that  the  "hangar  queen"  test  stand  can  repair  component 
i,  given  k  shortages,  is  then 


Pj(t.k) 


pkxT(t)ke-XT(t) 

k! 


(3) 


-(l-Pl)X(t) 


p^XT(t)VPiX  (t> 

— 


The  probability  that  the  hangar  queen  test  stand  can  repair  com 
ponent  i  at  time  t  is 


T 

P}(t) 


Et  — (1— P . ) X  (t) 

pj(t,k)  =  e 


X)p^T(t)k,-vT<t) 


k-0 


k-0 


(4) 


„  e-d-Pi)AT(t) 
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This  probability  is  used  in  the  simulation  to  determine  availability  of 
the  test  equipment  for  testing  certain  types  of  components.  Tests  are 
first  scheduled  on  the  hangar  queen,  utilizing  what  capability  exists, 

'  are  then  scheduled  on  the  remaining  pieces. 

SPARE  AIRCRAFT  COMPONENTS  TO  OVERCOME  QUEUES 

Additional  spare-parts  inventory  levels  are  determined  for  cases  in 
which  there  is  insufficient  service  capacity.  These  spare  parts  allow 
the  aircraft  to  remain  at  the  desired  level  of  mission  capability  even 
though  certain  components  are  not  repaired.  Clearly,  this  is  an  option 
only  for  a  limited  time  period,  since  it  is  never  possible  to  provide 
spares  for  insufficient  service  capacity  over  an  unlimited  time  span. 

The  approach  taken  is  to  determine  whether  the  average  service  capacity 
will  be  exceeded  over  the  time  horizon  of  interest,  and  when  it  is,  to 
"buy  out"  enough  low-cost  components  to  allow  the  average  service  capac¬ 
ity  to  match  the  average  service  demands.  Thus,  the  procedure  is  as 
follows : 

1.  Determine  the  average  service  capacity  of  each  server  type  for 
the  period. 

2.  Determine  the  average  service  demands  for  each  server  type  for 
the  period. 

3.  Determine,  for  those  server  types  with  insufficient  capacity, 
the  component  with  the  lowest  ratio  of  cost  to  test-time.  Buy 
spares  for  that  component  until  either  service  capacity  is  ade¬ 
quate  (assuming  that  each  spare  represents  a  test  that  does  not 
have  to  be  performed  during  the  period)  or  until  the  average 
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demands  for  the  component  for  the  period  are  covered.  Continue 
until  adequate  service  capacity  exists  for  each  server  type. 

These  "server  inadequacy"  spares  are  used  along  with  input  supply 
levels  as  the  starting  basis  for  supply-level  determination  in  the  pre¬ 
vious  sections. 


1 

-J  -  •-  -  -  -  * 
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Appendix  A 

CONVERGENCE  OF  THE  MARGINAL  ANALYSIS  APPROACH 

This  appendix  describes  conditions  of  convergence  for  the  basic  sup- 

ply  optimization  algorithm  of  Sec.  VII,  problem  (1). 

That  problem  is 

restated  here  as  follows: 

1 

minimize  CS  =  ^  S^ 

i  components 

1 

subject  to  P(Kjj  , S)  >  l 

(1) 

T 

S  >  S.  for  each  i 

i  -  io 

integer  . 

The  marginal  analysis  algorithm  attempts  to  find  an  " 

undominated 

'  solu- 

tion  to  this  problem. 

m 

An  undominated  solution  is  defined  by  Fox  (1966) 

as  one  for 

which 

m 

CS  <  CS*  =>  PCK^ ,S)  <  PCK^S*) 

(2) 

and 

'  m 

i 

CS  <  CS*  ->  PC^.S)  <  PO^.S*) 

(3) 

'  m 

.  J 
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That  is,  since  by  definition  PfK^.S*)  2:  a,  any  solutions  that  are  better 

than  S*  must  have  a  value  for  PfK^jS)  that  lies  between  P(K^,,S*)  and  a. 

If  PCK^S*)  =  a,  then  S*  is  an  optimal  solution.  Since  the  value  of  S* 

is  arrived  at  by  increasing  some  component's  supply  level  by  one  unit  at 

the  last  step,  we  know  that  P (K^ , S*)  is  likely  to  be  close  to  a. 

k-1 

Let  S*  be  arrived  at  on  the  kth  iteration  and  S  be  the  previous 
value  of  S.  Let  S  be  the  optimal  solution  to  problem  (1).  Then  we  also 
know  that 


P(KN,Sk  l)  <  a  <  P(Kn,S)  <  P(Kn,S*)  . 


(4) 


The  reason  that  logarithms  are  used  is  to  guarantee  convergence  to  an 
*  undominated  solution.  Marginal  analysis  will  provide  an  optimal  solu 

tion  to  the  problem 

minimize 


'  « 


subject  to  <)>(S)  >  ot 


S  >  0  for  all  i 


if  <j>(S)  is  separable  into  <KS)  =]C!tl1(s1)  and  each  of  these  is  a 
concave,  increasing  function.  See  Fox  (19  66).  Clearly,  the  function 
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P(KN,S)  *  II  pi(Q1*KN) 


(6) 


is  not  separable.  However,  we  note  that 


P(Kn,S)  >  a  <=>  An  PC^.S)  >  An  a  (7) 

and 

1 

PO^.S)  <  a  <=>  An  P(KN>S)  <  An  a  (8) 


t 


« 


when  P(lO,S)  is  an  increasing  monotonic  function  of  S;  therefore  an 
undominated  solution  to 


minimize 


E 


ciSi 


subject  to: 


An  P(K.t,S)  >  An  a 
N  ~ 

S  >  0 


(9) 


J 


is  an  undominated  solution  to  (1).  Furthermore, 

An  P(Kn,S)  -  An  IlPi(Q  K  )  -  E  An  P1(Q.K)  ,  (10) 

£  A  ^  1  w 

so  that  we  have  a  separable  problem.  It  remains  to  show  when 
An  P  (Q^)  is  concave.  Concavity  of  a  function  f(S1>  on  integer 
values  of  S±  is  defined  by 

f(s1)  -  f(S1  -  1)  >  f(S1  +  1)  -  f(Si)  .  (11) 

Let 

P  E  P  (QfKjj)  to  8^ow  dependence  on  . 

If  P  is  concave,  then  An  P  (Q^K^,S^)  will  also  be  concave. 

Now, 


Then 


,  QiKN«i 

p  <QilSi'si>  -  £  • 

k-0 


(12) 


pl«iKs'si>  -  plW1«B.s1  -  1)  -  PjCQi^.Sj)  , 


(13) 


1 
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and 


p  +  -  P  (QiVV  =  Pi<QiVSi  +  l) 


(14) 


Thus,  we  have  concavity  if 


Pi«iVSi  +  »  <  PiWlItN'Sl) 


(15) 


The  recursive  equations  for  the  two  distributions  in  Dyna-METRIC  are: 


pk+1  ,  MO  pk 
i  (k+1) 


ak+1  .  •  £  pk 

1  k+1  q  1 


Poisson 


Negative  Binomial 


(16) 


(17) 


where 


r  -  hill 

q-1  * 

q  =  variance-to-mean  ratio,  and 
p  =  q  -  1. 

For  each  of  these  distributions  we  see  that  if  k  >  X(t),  then  P 


,k+l 


<  P^  and  the  conditions  for  concavity  are  satisfied.  Thus,  the  separ¬ 
able  terms  of  the  constraint 


in  P  (Q^-V 
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are  concave  for  S.j+Q^K^  >  A^(t).  For  a  reasonable  value  of  confidence 
(that  is,  for  a  reasonable  value  of  o)  and  mission  non-capability  (K^) , 
we  can  therefore  expect  the  functions  to  be  concave  and  the  conditions 
for  convergence  to  be  satisfied. 
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Appendix  B 

THE  INDENTURED  SUBPROBLEMS 


This  appendix  describes  the  derivation  and  form  of  the  indentured 
subproblems  for  optimization  problem  (18)  in  Sec.  VII. 

Given  primary  assembly  stock,  we  wish  to  minimize 


The  function  B(Ssl)  represents  the  average  shortage  allocated  to  com¬ 
ponent  i  due  to  shortages  in  its  subcomponents.  As  shown  in  Sec.  V,  the 
worst  case  is  to  assume  that  the  sum  of  backorders  of  the  subcomponents 
contributes  directly  to  additional  primary  components  in  the  pipeline. 
(This  would  not  be  the  case  if  cannibalization  of  subcomponents  were 
performed.)  Again,  we  have  used  this  conservative  approach  because  gen¬ 
erally  the  much  lower  cost  of  subassemblies  allows  a  little  more  leeway 
in  "overbuying."  Under  this  assumption  we  have 

N  subassemblies 
B(S81,  -  £»,«,•>. 

j-1 


(3) 
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To  create  a  separable  problem  for  the  subassemblies,  we  use  a 


i  „siv 


linear  approximation  of  (Q^K^, S'1", Sax)  with  respect  to  B(SS1).  Thus 


si 


pl(<^i^N,Si  ,sSl) 


N  3P  (Q^.S^S  ) 

J=1  ^Bj(S j) 


i  -Si 

+  ?  (Q^.S^S  )  . 


(Bj(Sj) 


Vj” 


(4) 


Leaving  out  the  constant  terms  and  substituting  the  approximation  into 
problem  (4)  gives  the  revised  problem: 


min 


£cV  -  irZ- - 1 - jrE 

j  3  J  i  P  (Q1KN,Si,S81)  j 


i  -8i 

dPX(QK,S  S  ) 

- -  *  Bi(si} 

SW 


(5) 


Sj  >  0  for  all  j  . 


This  separates  into  subproblems  for  each  subcomponent  in  the  form 


a  8  8 

min  cjSj 


-  Y 


si 


V-T. 


3P  (Q1Kn,S1,S  ) 


,Si, 


i  P  (Q1KN,s1,s"*)  3B  (S") 
Sj  >  0  for  all  j  . 


•Vs°> 


(6) 


"1 


Let  (5j  be  the  value  of  the  term  in  brackets.  Then  we  have  the  problem, 


•in  c!sr  ’'Wsj) 


Sj>  °. 


which  is  easily  solved  (once  a  and  are  known)  by  increasing  from  0 
until  the  objective  function  no  longer  decreases.  The  result  is  the 
amount  of  subassembly  stock  and  completes  step  4.  The  value  of  (3^  must 
be  determined  through  perturbation  of  S_.  in  P  CQ^K^S^.S"’  )  except  for 
the  case  of  the  Poisson  distribution.  In  this  case  it  is  easy  to  show 


3P1(Q1KN,Si,S  ) 

*V>J> 


(X4(t)  +  B(S~)) 

(QiKN  +  V  * 


si  .AV8! 


j-(Xi(t)  +  B(S81))  •  (8) 


Furthermore,  for  the  Poisson  case,  the  solution  to  (7)  can  be  converted 

g 

to  a  ready-rate  criterion  on  S  ^  .  The  solution  to  (7)  requires  that 


c j(S *  +  1)  -  CjS®  -  y  ej[Bj(S®+  1)  -  Bj(S®)]  >  0  . 


For  the  Poisson  distribution,  the  difference  in  brackets  is  merely 
[RRj (S8)  -  1],  This  gives  us  the  criterion 


YSjtRRj  (S®)  ~  1]  <  c® 


(10) 


which  implies  that  the  solution  to  (7)  can  be  obtained  by  increasing  the 
stock  of  component  j  until  it  has  a  ready  rate  satisfying  this  cri¬ 
terion. 
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